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The Diffraction of X-Rays by Liquid Benzene-Cyclohexane Mixtures 


A. EISENSTEIN AND N. S. GINGRICH 
University of Missouri, Columbia, Missouri 


(Received February 2, 1942) 


The x-ray diffraction pattern of a liquid mixture of benzene and cyclohexane is investigated, 
using both the photographic method and the Geiger-Mueller tube counter method. The fine 
structure on the main peak, and a secondary peak as reported by Bell and Davey, cannot be 
obtained. It is concluded that these features do not exist in this pattern. 





INTRODUCTION 


N a recent article, Bell and Davey! reported 

that they found some very unusual fine 
structure in the main diffraction peak from 
benzene-cyclohexane mixtures. A repetition of 
this work is highly desirable for three reasons: 
first, no liquid x-ray diffraction pattern has 
heretofore shown fine structure of this type; 
second, this particular mixture has been worked 
upon by Ward,? by Murray and Warren,’ and by 
Bell and Davey,' and all three reports of these 
workers are in disagreement; third, the extent of 
succeeding arguments based on the results of 
Bell and Davey is so great that it is desirable to 
examine critically the experimental basis for 
these arguments. The present authors have 
made an experimental study of the diffraction 
of x-rays by benzene, by cyclohexane, and by 
one mixture of these liquids, using both the 
photographic and the Geiger-Mueller counter 
methods of detection of the x-rays. Particular 
attention was given to the controversial matter 
of fine structure in the main peaks of the mixture. 





*P. H. Bell and W. P. Davey, J. Chem. Phys. 9, 441 


(1941). 
*H. K. Ward, J. Chem. Phys. 2, 153 (1934). 


*G. E. Murray and B. E. Warren, J. Chem. Phys. 7, 
141 (1939), 
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EXPERIMENTAL 


Benzene and cyclohexane of high purity 
were purchased, and no further purification of 
the liquids was attempted. The index of refrac- 
tion of the benzene used here was measured and 
found to be 1.5008 at 20°C while that used by 
Bell and Davey was given as 1.5009. The index 
of refraction of the cyclohexane used here was 
measured and found to be 1.4263 at 20°C, and 
that used by Bell and Davey was given as 1.4263. 
The mixture studied was approximately 50 per- 
cent benzene—50 percent cyclohexane by vol- 
ume, and its index of refraction was 1.4597. 
The cell holder consisted of a brass cylinder 1.1 
cm long and 1.9 cm diameter with thin Cello- 
phane windows clamped on the ends. After the 
various liquids were contained in the cell for 
several days, tests on the index of refraction 
were made. In the cases of benzene and of cyclo- 
hexane alone, no difference in the index of refrac- 
tion was ever obtained. In the case of the 
mixture, no difference greater than 0.0001 was 
obtained when adequate care in rinsing the cell 
was taken and when precautions were taken 
against evaporation through the filling tube. 

Two cameras were used in the photographic 
studies here reported, a small camera of 9.20-cm 
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Fic. 1. X-ray diffraction patterns of a 50-50 mixture of liquid benzene and 
cyclohexane, showing intensity as a function of distance on the films in cm. One 
cm corresponds to about 1° 16’. Curve A taken with Mo Ka radiation; curve B 
with Mo radiation filtered through 0.08 mm Zr equivalent; curve C with Mo radi- 
ation filtered through the above Zr and 0.132 cm Al. 


radius and a large camera of approximately 44.6- 
cm radius. Using the small camera, no indication 
of fine structure on the main peak could be 
obtained with either pure liquid or with the 
mixture, but it was observed that in pictures 
taken with filtered radiation a weak peak ap- 
peared at nearly half the angle for the main peak 
whereas with crystal reflected Mo Ka radiation, 
no secondary peak could be observed. Further 
work on the mixture was done using the large 
camera and three different types of radiation, 
first, crystal reflected Mo Ka, second Zr filtered 
Mo radiation, and third, Zr and Al filtered Mo 


radiation. Intensity curves obtained from micro- 


photometering typical films are shown in Fig. 1. 
The intensity, computed in the usual manner, is 
plotted as a function of distance on the film from 
the center of the film. In the lower curve, taken 
with Mo Ka crystal reflected radiation, it is to 
be noted, first, that no fine structure appears on 
the main peak, and second that no secondary 
peak at roughly half the distance (or angle) of 
the main peak can be detected. This is in con- 
tradiction to the results of Bell and Davey who 
state that in all cases a small intensity peak was 
observed at about 4°, and which had an intensity 
approximately one-eighth of the high intensity 
peaks. 
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Fic. 2. A microphotometer trace for the benzene-cyclohexane mixture, using Mo Ka 
radiation, showing the main diffraction peak on each side of the direct x-ray beam. The 
three arrows indicate the approximate position of the three fine structure peaks as re- 


ported by Bell and Davey. 


The total horizontal angular divergence of the 
incident beam was measured photographically to 
be slightly less than 7 minutes, a value which 
compares favorably with that of Bell and Davey. 
The density of the film at the peak was of such 
value that the fine structure as reported could 
easily have been detected. Furthermore the 
region of the film over which the fine structure 
should have appeared is more than half a cen- 
timeter, the general location of the three peaks 
being indicated by the three spurs in Fig. 1. 
Both the film and the microphotometer are 
easily capable of showing up structure of this 
type. Figure 2 shows a portion of a micro- 
photometer trace typical of the one from which 
the lower intensity curve of Fig. 1 was obtained. 
The gradual increase of intensity at small angles 
is due to air scattering, and the cut-off at still 
smaller angles is due to a lead strip to suppress 
the direct beam. 

In an attempt to find a possible explanation 
for the secondary peak at about half the angle 
of the main peak, two other patterns were 
obtained. These are shown in Fig. 1, the middle 
curve being taken with Mo radiation filtered 
with ZrO, of 0.08 mm Zr equivalent thickness, 
the upper curve with Mo radiation filtered with 
the same Zr filter and 0.132 cm of Al. In both 
cases, the total horizontal angular divergence of 
the incident beam was measured photographi- 
cally to be 10 minutes, and the same camera was 


used for all three curves of Fig. 1. The origin of 
the extra peak at nearly half the angle of the 
main peak is unquestionably that which is given 
in detail elsewhere.* Residual continuous radi- 
ation forming the spurious peak in the Zr filtered 
pattern is greatly enhanced in the Zr+Al 
filtered pattern due to the greater absorption of 
0.710A as compared to that of the continuous 
radiation at about 0.4A. This rough effective 
wave-length is that obtained with the x-ray tube 
operated at about 40 kv peak, whereas Bell and 
Davey, operating their tube at 59 kv peak, 
obtained the secondary peak at an angle smaller 
than that obtained in the present work. This is 
qualitatively in accord with the above possible 
explanation of the origin of their secondary peak. 
In the use of the balanced filter method, extreme 
care must be taken to be sure that a perfect 
balance is obtained over the entire region 
studied, especially when compounds of the 
balancing elements are used. 

Using a Geiger-Mueller tube counter, an am- 
plifier, a scale of sixteen and a Cenco counter, a 
further study was made of the diffraction of 
x-rays by the same liquid benzene-cyclohexane 
mixture as was used above. Rocksalt reflected 
Mo Ka radiation was passed through two slits 
as shown in Fig. 3. The total horizontal angular 
divergence of the incident beam, determined 


4R. Q. Gregg and N. S. Gingrich, Rev. Sci. Inst. 11, 
305 (1940). 
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Fic. 3. Geometrical arrangement using the Geiger-Mueller counter method. 


largely by the crystal, was measured photo- 
' graphically and found to be 10 minutes. The 
maximum possible divergence of the diffracted 
beam through the last two slits is 26 minutes or 
about the same as that of Bell and Davey’s 
Soller slit system. Many runs were taken, par- 
ticularly over the top of the main peak and the 
average of these runs is shown in Fig. 4, together 
with a reproduction of one of Bell and Davey’s 
curves at roughly the same temperature and 
mixture. The maximum counting rates in the 
two cases are roughly 120 per minute in the 
previous work, and 470 per minute in the present 
work. The total number of counts per point in 
the present work is from 12,000 to 16,000. The 
probable error in the determination of each 
point was calculated according to the expression 


p.e. = +0.67N/(NT)} 


where JN is the rate of counting, and NT the 
total counts. This probable error is indicated on 
the lower curve of Fig. 4. It is seen that though 
the points do not fall on a perfectly smooth 
curve, the curve drawn does pass through some 
portion of the limits of the probable errors. In 
no case could the results here obtained show a 
fine structure anywhere near as large as that 
reported by Bell and Davey. In preliminary runs 
in which only a few hundred counts per point 
were taken, irregularities were obtained which, 
if no further work was done, might conceivably 
have been interpreted as fine structure, though 
never, in the present work, as great as that in 
question. Furthermore, these irregularities ap- 
peared quite at random. It was found, for 
example, that reliable curves could be taken only 
after midnight and before 7 a.m., during which 


time the line voltage was most steady. Even 
assuming ideal experimental conditions, however, 
the results, using a Geiger counter, can only be 
statistical in nature, and fluctuations always 
appear more prominent on a flat or nearly flat 
portion of the curve. 


CONCLUSIONS 


In a photographic method of recording the 
x-ray diffraction pattern of a liquid benzenc- 
cyclohexane mixture, especially designed to 
show any possible fine structure, no significant 
trace of fine structure could be observed such as 
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Fic. 4. Lower curve shows counts per minute as a 
function of scattering angle, 26, as determined in the 
present work. The probable error is indicated by vertical 
arrows. Upper curve is one reproduced from Bell and 
Davey’s work. 
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was reported by Bell and Davey. Further studies 
were made, demonstrating the fact that impure 
radiation can give rise to a secondary peak at 
small angles. It is suggested that it may be 
possible that this is the origin of Bell and 
Davey’s secondary peak, caused by a slight 
unbalance of their balanced filters. 

Using the Geiger-Mueller counter method, as 
Bell and Davey did, but with pure radiation and 
with geometrical conditions at least as good as 
in their work, no fine structure could be obtained. 
Slight fluctuations of the observed points about 
a smooth curve are not large enough to establish 


fine structure, and these fluctuations are con- 
sidered to be wholly statistical. 

It is concluded, further, that the extensive 
deductions made by Bell and Davey are without 
experimental foundation. 
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The spectrum of powdered AgCN has been studied in the region 2500-1800 cm™. Intense 


absorption at 2178 cm™ is attributed to an oscillation involving the CN group. The CN fre- 


NVESTIGATIONS of the infra-red and 

Raman spectra of the organic cyanides have 
revealed the existence of a vibrational frequency 
of approximately 2245 cm! which is charac- 
teristic of the CN groups, while similar studies 
of the organic isocyanides indicate that the NC 
group has a characteristic frequency of 2146- 
2161 cm-'. The simpler inorganic cyanides, 
NaCN and KCN, have CN frequencies of 
approximately 2080 cm-', while in the complex 
cyanides such as KAg(CN)2 and K;Fe(CN)6. 
-3H.O the characteristic frequencies are in the 
neighborhood of 2150 cm-'.4? It has been sug- 
gested by several workers* that in inorganic 
cyanides the CN frequency increases when there 





‘Summary of the Raman data: K. W. F. Kohlrausch, Der 
Smekal-Raman Effekt (J. Springer, Berlin, 1931), p. 320; J. 
H. Hibben, The Raman Effect and Its Chemical A pplications 
(Reinhold, New York, 1939), p. 452. 

Infra-red data: W. Gordy and D. Williams, J. Chem. 
Phys. 3, 664 (1935); 4, 85 (1936); F. K. Bell, J. Am. Chem. 
Soc. 57, 1023 (1935). 

* P. Krishnamurti, Ind. J. Phys. 5, 651 (1930); R. Samuel 
and M. J. Khan, Zeits. f. Physik 84, 87 (1933); Volkensh- 
tein, Acta Physicochim. (U.S. S. R.) 7, 315 (1937). 


quency for AgCN is almost 100 cm™ higher than the CN frequencies of NaCN and KCN. 


is an increase in the homopolar nature of the 
M-CN binding. 

It seemed desirable to test this hypothesis by 
determining the CN frequency for a simple 
cyanide in which the M—CN bond possesses 
some homopolar properties. As the existing data 
provide only comparisons of NaCN and KCN 
with double cyanides and with complex cyanides 
of the type mentioned above, the study of the 
infra-red spectrum of AgCN for comparison with 
NaCN immediately suggested itself, since the 
homopolar properties of the silver salts are well 
known.! 

The silver cyanide used in the present work 
was prepared by adding an equimolecular amount 
of potassium cyanide solution to an approxi- 
mately 0.1 solution of silver nitrate with 
constant and vigorous stirring. After precipita- 
tion was complete, the supernatant liquid was 
decanted and the precipitate was placed in a two 


4L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, 1940), p. 73. 
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and a half foot large glass test tube. A fine jet of 
water was then directed into the precipitate. In 
this manner the precipitate was continuously 
washed for two days. The precipitated silver 
cyanide was then collected and dried. 

The spectrometer used was a rocksalt prism 
instrument of the Littrow type. The samples 
used in making the absorption measurements 
consisted of layers of powder deposited on 
fluorite plates. The thickness of the absorbing 
layers varied from 0.02 to 0.05 mm in different 
samples. As the mean diameter of the particles 
was approximately three microns, considerable 
scattering took place in the spectral region 
studied. In order to minimize this scattering the 
powder was immersed in carbon tetrachloride in 
some of the cells.5 

Figure 1 illustrates the type of transmission 
curve obtained with powdered AgCN. Curve A 
represents the transmission of a typical sample, 
the dotted curve giving the estimated back- 
ground due to scattering. Curve B shows the 
transmission of AgCN when corrected for scat- 
tering. It will be noted that there is an intense 
absorption band with maximum at 4.59u. This 
band was observed for a number of different 


TABLE I. The CN oscillational frequencies in 
various compounds. 











Compound Infra-red Raman 
NaCN 2080 GW 2085 K 
KCN 2080 GW 2081 P-SG 
AgCN 2178 PS -_ 
Ba(CN)2 2080 GW —- 
Hg(CN). wos 2192 K 
CH;NC 2183 GW 2161 B 
C:H;sNC 2159 GW 2146 B 
CH;CN 2285 GW 2250 P-SG 
C:H;CN 2258 GW 2246 DK 








K—Krishnamurti, Ind. J. Phys. 5, 633 (1930). 

P-SG—Pal and Sen Gupta, Ind. J. Phys. 5, 13 (1930). 

DK—Dadieu and Kohlrausch, Wien Ber. 139, 165 (1930). 

B—Bhagavantam, Ind. J. Phys. 5, 48 (1930). 

GW—W. Gordy and D. Williams, J. Chem. Phys. 3, 664 (1935); 
4, 85 (1936). 

PS—Present study. 


5 Carbon tetrachloride has no absorption bands in the 
range 2500-1800 cm which was studied. 
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Fic. 1. The transmission of silver cyanide in the 2500- 
1800 cm“ region. 





cells, and the frequency of maximum absorption 
was found to be 2178 cm“. 

Table I summarizes the pertinent infra-red 
and Raman data for several types of cyanides. 
It will be noted that the CN frequency in AgCN 
is almost 100 cm higher than for NaCN and 
KCN. This result is in accord with the sug- 
gestion that increasing homopolarity is accom- 
panied by increases in the CN frequency. It 
might be mentioned in passing that the AgCN 
frequency actually falls in the range of the 
organic isocyanide, NC, frequencies, although 
many factors render unjustifiable any suggestion 
of similarity in structure from spectroscopic 
evidence. Included in the table are data for two 
divalent cyanides, Ba(CN)2 and Hg(CN):. In 
this case also the compound having the higher 
degree of homopolarity has the higher CN 
frequency. It will be observed from the table 
that none of the inorganic cyanides yet studied 
has as great a CN frequency as that characteristic 
of the organic cyanides. 
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The second virial coefficients of ethane, propane, 
n-butane, n-heptane, ammonia, methyl chloride, and the 
freons are computed from available experimental data. 
The causes for sizeable errors in second virial coefficients 
are considered. At temperatures above the critical, the 
second virial agrees with the theorem of corresponding 
states. Below the critical temperature, molecules with 
dipoles have unusually large virials and the values of their 
reduced dipole moment, »/(T-V-)', determine the dis- 
crepancy. The data for isomeric hydrocarbons show that 
the second virial is not sensitive to the shape of the mole- 
cule. This makes it impossible to determine the exact laws 
of intermolecular interaction from the temperature vari- 


ation of the second virial. The imperfections of a gas are 
considered to arise from the presence of double molecules 
which exist for the duration of a collision. The equilibrium 
constant for the formation of double molecules is related 
to the second virial and its temperature variation gives 
the entropy AS and the energy AE of their formation. 
These may be interpreted qualitatively in terms of the in- 
termolecular forces. The second virial coefficient for all sub- 
stances fits the equation: B(T)=6b,(1—c exp (—AE/RT)) 
where ba, c, and AE are taken as constants. Here 
bac=exp (AS/R). A corresponding states equation is given 
for estimating second virials when no experimental data 
are available. 





N this paper we investigate the significance of 

second virial coefficients for complex mole- 
cules. The treatment is divided into three parts: 
(1) the determination and probable accuracy of 
experimental values; (2) the qualitative inter- 
pretation of second virials from the theorem of 
corresponding states; and (3) a quantitative 
interpretation in terms of the energy and entropy 
of formation of instantaneous double molecules. 
Unfortunately it is difficult to distinguish the 
characteristic properties of the molecules because 
the second virials may have errors of ten cc/mole. 
At temperatures above the critical, the second 
virial agrees with the theorem of corresponding 
states. Below the critical temperature, molecules 
with dipoles have unusually large virials and the 
value of their reduced dipole moment, u/(T.V-)}, 
determines the discrepancy. An examination of 
Roper’s second virial coefficients for the isomeric 
hydrocarbons shows that the second virial 
coefficient is not sensitive to the shape of the 
molecule. 

For the practical interpretation of the second 
virial coefficients of complex molecules, it is 
useful to exploit an idea of Fowler and Guggen- 
heim.2 Consider an imperfect gas as an equi- 
librium mixture of single, double, triple, etc. 
molecules. The double, triple, etc. species are 





'E. E. Roper, J. Phys. Chem. 44, 835 (1940). 
* R. H. Fowler and E. A. Guggenheim, Statistical Thermo- 
dynamics (Cambridge University Press, 1939), Chap. 7. 


complexes which exist only for the duration of a 
collision, but their number is determined by the 
usual type of chemical equilibrium constant. 
The equilibrium constant for the formation of 
double molecules is related to the second virial 
coefficient, and its temperature dependence gives 
the energy and entropy of their formation. This 
energy and entropy of formation of double 
molecules may be roughly interpreted in terms 
of intermolecular forces. To determine these 
forces quantitatively, we would require more 
accurate experimental data. 


I. THE DETERMINATION OF SECOND 
VIRIAL COEFFICIENTS 


The equation of state of a gas at moderate 
pressures can be written in the form: 


pV=RT+pB(T)+P°C(T)+---. (1) 


This is called the virial form of the equation of 
state since Clausius first derived it by equating 
the total kinetic energy of the molecules to the 
total virial of the system, —})0; (x;F2:+yiF uv: 
+2;F2z;). Here x; is the x coordinate of the ith 
particle, and Fz; is the x component of the force 
acting on the ith molecule, etc. The term pV 
arises from the forces exerted on the walls of the 
containing vessel, RT comes from the kinetic 
energy of the molecules, and all of the other 
terms arise from the mutual interactions between 
the molecules. B(T) is called the second virial 
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Fic. 1. The four types of curves obtained by plotting 
the experimental values of Q= V—RT/p as a function of 
pressure. 


coefficient, C(T) is the third virial coefficient, 
etc. The second virial coefficient is particularly 
useful for our purposes because it arises only 
from collisions between two molecules at a time 
and is related rigorously to the energy of inter- 
action between pairs of molecules. 

The easiest way of determining B(T) from ex- 
perimental p, V, T dataistoplotQ=(pV—RT)/p 
as a function of pressure. Since 


Q=(pV—RT)/p=B(T)+C(T)pt+---, (2) 


the second virial coefficient is the limit of Q as 
the pressure approaches zero: 


a) 
> ‘ 


The third virial coefficient is the slope of Q at 
zero pressure, etc. When we plot actual experi- 
mental data in this manner, we get four types of 
curves as shown in Fig. 1. Curve (1) is typical 
of substances above their critical temperature 
where the third and higher virial coefficients do 
not become important until rather high pressures 
are reached. Curve (2) is typical of substances 


B(T) =lim (3) 
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below their critical temperature, the sharp break 
in the curve occurring when the pressure 
approaches the vapor pressure. Curve (3) could 
correspond to a gas which dissociates at low 
pressures to give m molecules. In such a case, 
PV approaches RT at low pressures and 
Q=(pV—RT)/p approaches infinity. However, 
curves of the form (3) can correspond to even 
small experimental error in the determination of 
the volume of the apparatus. Curve (4) must 
always correspond to errors of this nature. 
When the experimental details are available, 
it is sometimes possible to make a single cor- 
rection to the volume of the apparatus such 
that the curves for Q now resemble either curves 
(1) or (2) for all temperatures. Otherwise it is 
necessary to extrapolate Q to zero pressure from 
the experimental data at moderate pressures 
where a small error in the volume of the appa- 
ratus has only a slight effect. Either of these 
processes naturally reduces the accuracy of the 
second virial coefficient. Unfortunately many of 
Sage and Lacey’s measurements of p, V, T for 
hydrocarbons lead to curves of the form of (4) 
and are therefore unsuitable for the determina- 
tion of second virial coefficients. This behavior 
was also manifested by the data for steam 
which Keenan and Keyes tabulated in the first 
edition of their “Steam Tables” but this has 
been corrected in their later editions.* We believe 
that accurate values of the second virial coefh- 
cients can best be obtained from specially 
constructed low pressure apparatus such as 
Roper! used in his excellent work on the olefins. 

The second virial coefficients which we com- 
puted from the original experimental data are 
given in Table I. We were surprised at the lack 
of accurate low pressure p, V, T measurements 
of the paraffins. Only Deschner and Brown’s' 
measurements for propane went to pressures as 
low as one atmosphere. The ethane, butane, 
and heptane results as well as the Beattie, Kay, 
and Kaminsky® values for propane were extra- 
polated from pressures as high as 25 atmospheres. 
The virial coefficients for steam and for ammonia 


3 J. H. Keenan and F. G. Keyes, Thermodynamic Prop- 
erties of Steam (Wiley, 1936). 

4W. W. Deschner and G. G. Brown, J. Ind. Eng. Chem. 
32, 836 (1940). 

5 Beattie, Kay, and Kaminsky, J. Am. Chem. Soc. 59, 
1589 (1937). 











(1S 
§ 
(19 








oC — = ~ 


im 
rst 
jas 
ve 


liy 
as 


ym- 
are 
ack 
nts 
n’s! 
as 
une, 
ayy 
tra- 
res. 
onia 


we 


Prop- 
‘hem. 


pe 59, 








should be quite accurate. The data for methyl 
chloride seem to be quite reliable. Of the freons, 
CCI,.F.—Freon 12, is the least accurate since the 
temperature measurements were uncertain to 1°. 
Second virial coefficients for other boiler fluids 
and refrigerants can be obtained from a study 
of similar commercial publications. The virials 
for the olefins are given directly in a paper of 
Roper! and are therefore not reproduced here. 
It is difficult to obtain the second virial 
coefficients to within 10 cc/mole. For this 
accuracy we should know Q(7T)=V—RT/p to 
within 10 cc/mole. At room temperature, V is 
of the order of 24,000 cc/mole when is of the 
order of one atmosphere. Thus the density of 
the gas, the pressure in the chamber, and the 
absolute temperature must all be measured to 
within one part in 2400. This is feasible with 
carefully built apparatus. Often we lose sight 
of the fact that the molecular weight and the 
gas constant must be known to the same accu- 
racy, and the difference between two published 
sets of virial coefficients may be due to the 
different values of these constants which were 
used. Thus second virial coefficients for organic 
molecules published prior to 1937 when the 
atomic weight of carbon was thought to be 12.00 
may differ very considerably from those in 
which the present atomic weight of carbon, 
12.01 is used. However, for the saturated 
hydrocarbons where the older molecular weights 
were used,5-8 the errors due to this cause would 
be inappreciable since all of the measurements 
were at high pressures. When Amagat units 
are used, molecular weights are not required and 
it is only necessary to know the gas constant R 
with reasonable accuracy. In recent years there 
has been some controversy as to whether the 
ice point is 273.16°K or 273.13°K. In the former 
case, R=82.054 atmos. cc per mole per deg. 
instead of 82.060 and the difference in B(T) is 
approximately 3 cc/mole. Holborn and Otto 
(1925) and Kammerlingh-Onnes (1910) avoided 
serious errors in the virials for argon by using 
Amagat units and not relying on the faulty 





* Beattie, Hadlock and Poffenberger, J. Chem. Phys. 
3, 93 (1935). 

’ Beattie, Simard, and Su, J. Am. Chem. Soc. 61, 26 
(1939), 


i937 Beattie, and Kay, J. Am. Chem. Soc. 59, 1587 
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atomic weights available then. Often two good 
sets of data such as Deschner and Brown’s and 
Beattie, Kay, and Kaminsky’s for propane 
(see Table I) differ by a constant independent 
of temperature. Presumably this is because of 
small differences in the constants used or in the 
calibration of the thermometers. Joule-Thomson 
measurements are inherently more accurate than 
second virial coefficients because they do not 
require taking the difference of two large 
quantities. Also, Joule-Thomson coefficients are 
just as easy to interpret as the virials. 

Often second virial coefficients are estimated 
from empirical equations of state such as the 
Beattie-Bridgman. This method has the ad- 
vantage of being quick and easy, but it is 
sometimes unreliable. The constants for empirical 


TABLE [. 








Second virial coefficient, B(T), in cc/mole 





Propane 


T°K Ethane re BKK® DBe mF n-Butaneé 
298.16 —191 303.16 —384 423.1 —324 
323.16 —160 348.16 —293 448. is —288 
348.16 —139 369.97 —260 473.16 —256 
373.16 —122 373.16 —247 —256 498.16 —224 
398.16 —108 398.16 —211 —224 523.16 —198 
423.16 — 94 423.16 —183 —197 548.16 —176 
448.16 — 86 448.16 —160 573.16 —159 
473.16 —- 77 473.16 —139 —155 
498.16 — 71 498.16 -121 
523.16 — 60 523.16 —109 
526.38 —108 
548.16 — 96 
570.46 — 89 
Methyl 
Chlorides 
T°K n-Heptane* T°K Ammonia‘ T°K CH;Cl 
548.16 —569 273.16 —345 238.71 —764 
573.16 —503 298.16 —261 249.82 —668 
598.16 —449 323.16 —209 255.38 —637 
623.16 —399 373.16 —142 283.16 —500 
423.16 —101 310.94 —401 
473.16 — 75 338.71 —320 
523.16 — 58 366.49 —265 
573.16 — 44 394.27 —214 
422.05 —184 
449.82 —155 
Freon—114 Freon—21* Freon—12/ Freon —113* 
T°K CCI:sF CHChF CChF+ CChF —CCIF2 
238.71 —1150 —766 —637 
244.27 —604 
249.82 —1070 —734 —586 
255.38 —560 
283.16 — 862 —616 —454 — 1080 
310.94 — 720 —528 —347 — 944 
338.71 — 617 —446 — 795 
366.49 — 523 —403 — 737 
394.27 — 459 —354 — 654 
422.05 — 381 —310 — 561 
449.82 — 328 —271 — 495 








@ See reference 6. 

+ See reference 5. 

© See reference 4. 

4 See reference 7. 

¢ See reference 8. 

4 C. H. Meyers and R. S. Jessup, Refrig. Eng. 11, 345 (1925). 

9 Methyl Chloride (E. I. du Pont de Nemours and Co., The R. and 
H. Chemicals Department, 1940), fifth edition. 

4 Thermodynamic Properties of Trichloromonofluoromethane (Kinetic 
Chemicals, Inc., 1938). 

‘ Thermodynamic Properties of Dichloromonofluoromethane (Kinetic 
Chemicals, Inc., ). 

j Thermodynamic Properties of Dichlorodifluoromethane, Circular No. 
12, Am. Soc. Refrig. Eng. (1931 

k Thermodynamic Properties "of Trichlorotrifluoroethane (Kinetic 
Chemicals, ie. 1938). 











equations of state are usually made to fit high 
pressure measurements and therefore there is no 
guaranty that the equation really fits at low 
pressures. Furthermore, Michels and Michels® 
have shown from their work on CO, that there 
is an intrinsic difference between the virial 
coefficients which fit low pressure and those 
which fit high pressure experiments. This be- 
havior can be explained by the presence of 
exponential terms such as exp (—cV) in the true 
equation of state. Such a term is asymptotically 
zero for large volumes and therefore no trace of 
its existence would appear in the low pressure 
measurements. However, at high pressures, 
exp (—cV) can be expanded in an ordinary 
power series to give the anomalous virial coeffi- 
cients. On both theoretical and experimental 
grounds, such an exponential term might exist 
in the internal energy U and therefore appear 
in the equation of state. The exponential form 
of the quantum-mechanical exchange integrals 
together with the volume dependence of the 
internal energy of crystals suggests that: 


U=U)-a/V—d exp (—vV) 


and this form is in good agreement with the 
experimental observations of Hildebrand!’ and 
others. 

The usefulness of second virial coefficients for 
theoretical purposes depends on the rigorous 
statistical mechanical relation: 


N 
B(T)= 
167? 





f- . ft —exp CRE 


where 
dQ=sin 6déd¢dy 
0<0<2, 0<G@<27, 0<Y<27 
dr =dxdydz. 


Here E is the energy of a pair of molecules 
expressed as a function of the distance between 
the molecules and the relative orientations. In 
carrying out the integrations we can suppose 
that one molecule is held fixed at the center of 
a coordinate system and the second molecule is 
allowed to traverse all space with its center at 


9 A. Michels and C. Michels, Proc. Roy. Soc. A160, 348 
(1937). 
10 J. H. Hildebrand, Solubility (Reinhold, 1936), second 
edition. 
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x, y, 2. The angles 6, ¢, and y are the Eulerian 
angles which the second molecule makes with 
respect to the first. Since E does not involve the 
temperature, we may regard Eq. (5) as an 
integral equation for E as a function of the 
coordinates in terms of B as a function of 
temperature. The remainder of this paper is 
concerned with the solutions of this equation. 

The potential energy of interaction of spheri- 
cally symmetrical non-polar molecules can be 
approximated by the form: 


E(r) =4Enl(r0/r)”— (ro/r)®]. (4) 


Here E,, is the maximum energy of attraction of 
the two molecules and 7 is the collision diameter 
corresponding to E(r9) =0. Lennard-Jones carried 
out the integration of Eq. (5) with this form of 
the potential energy and obtained: 


B(T) = boF 12(Em/kT). (5) 


Here bp = 227Nyr,'/3, is the van der Waals constant 
for the collision diameter 79. The function 
Fi(E,/kT) is a power series which can be 
calculated readily for any desired value of 
En/kT. Table II shows the values of Fy. from 
which it is possible to graph the function and 
interpolate to obtain the values of Fj. for any 
other values of the argument. 

Recently, Stockmayer showed how to calculate 
the second virial coefficient for spherical mole- 








TABLE II.t 
Em/kT Fi2(Em/kT) = Em/kT Fi2(Em/kT) Em/kT — F12(Em/kT) 
0025 ~—-.3247 .1986 .2468 .6280 —1.0624 
0158 4949 .2079 .2236 .7046 —1.3398 
0315  .5258 .2500 1154 .7906 —1.6668 
0628 5148 .3147 —.0621 9953 —2.5171 
0791 4945 3296 —.1044 1.1025 —3.0064 
0995 4617 .3614 —.1969 1.3225 —4.1171 
1253. «.4128 3962 —.3012 1.5625 —5.5184 
1577 3433 4988 —.6237 2.2500 —11.0328 
Em/kT  Hi2(Em/kT) J12(Em/kT) K12(Em/kT) L12(Em/kT) 
1225 —.0072 —.00001 
.1555 —.0113 —,.00002 
.1936 —.0173 —.00004 
.2025 —.0189 —,.00005 
.2500 —.0288 —.00011 —.00 X10- —.000 X1075 
.3025 —.0424 —,.00023 —.01 X1074 —.001 X1075 
.3600 —.0607 —.00045 —.03 X1074 —.002 X1075 
.4900 —.1174 —,.00152 —.19 X1074 —.018 X1075 
.6400 —.2137 —.00453 —.88 X1074 —.135 X1075 
.7225 —.2844 —.00758 —1.82 X1074 —.348 X1075 
.8100 —.3731 —.01237 —3.64 X1074 —.845 X1075 
1.0000 —.6324 —.03160 —13.85 X10-4 —4,698 X1075 
Em/kT M12(Em/kT) N12(Em/kT) 
.4900 —.001 X 1076 
.6400 —.017 X10-6 —.02 X107 
.7225 —.052 X10-6 —.06 X1078 
.8100 —.151 X1076 —.22 X10 
1.0000 —1.195 X1076 —2.38 X1078 








] + Many of the values of Fi2(Em/kT) were given to us by Lennard- 
ones. 
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cules in which a dipole yu is imbedded. In this 
case he assumes that the potential energy 
between pairs of molecules has the form: 


E(r) -48-|(~) -() ] 


pe 
——[2 cos 6; cos 62—sin 6; sif @.cos¢]. (6) 
r3 


Here, of course, 7) and E,, have a slightly different 
interpretation from those in the equation for 
non-polar gases. Stockmayer™ makes most of 
his calculations on the assumption that either s 
is 24 or infinite. We believe that it is more 
reasonable to take s=12 which is in better 
accord with the repulsive forces in non-polar 
gases. Then, the second virial coefficient becomes: 


B(T) = bol Fi2(Em/RT) +x i2(Em/kT) 
+2x7S12(Em/RT) +2°K 12(En/RT) 
+2'L10(Em/k T) +x°Mi2(En/k T) 

+x°Ni2(En/kT)+-+-+]. (7) 


Here x=y'/E,2ro® is the significant variable 
which gives the effect of the dipole. We can 
express Hy, Ji2, Kio, Liz, etc. in terms of the 
derivatives of Fi2(E,,/kT). Remembering that 
extensive tables’ are available for 


Fy. 


d 
Gix(E,,/kT) = r—_(—), 
dT\ T 


we obtain: 


1/5 
Mh.= ~-(-FntGu) (En/kT), 


1 
Ji2= ——(3F12+2Gi2)(Emn/RT)', 
200 


29 5 
Ky,= — | (Fat Gu (Em/kRT)* 
88,200L \4 


+ (3Fi2+2Gi2)(Em/kT)® 
5 
= ({F'n+6's) (E./eT) (7c) 


4 W. H. Stockmayer, J. Chem. Phys. 9, 398, 863 (1941). 
” J. O. Hirschfelder, R. B. Ewell, and J. R. Roebuck, 
J. Chem. Phys. 6, 205 (1938). 
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11 3 
Liy.= — aa (27 +562) (En/kT)® 
793,800 2 


+(3Fit2G—2F 


3 
~56'x) (B=/aT)'| 


312 5 En\' 
My2.= — | (“Fn +2Gu) (—*) 
77(ALS)LA2 kT 


5 En\'* 
+ (1172+ 861-2262) (—) 
2 kT 


+ (12% +86..- 11 F’12—8G' 12 


5 } ° 
+Spruta’s)(=))] 
4 kT 
— 1.045,1097 713 E,.\° 
N= ————| (P+ 56) (=) 
12! 2 kT 


13 
oe ( 14F).+ eat ale 


E,, \ 
-56'»)(—~.) + (12%. 
kT 


+ 8G12—14F':2—10G6'12 


13 5 Em\" 
+7F'nt36"n)(——) |. 
4 2 kT 


The single and double primes indicate first 
and second derivatives with respect to E,,/kT. 
In Table II the numerical values of these 
functions are tabulated, so that it is a simple 
matter to obtain the effect of a dipole moment 
on the second virial coefficient. Fitting the 
second virial coefficient of water, we find x= 10, 
En/k=400°, and b)=22.9 cc. When T=400°K 
the non-polar contribution, boF\.(En,/kT), is 
only 17 percent of the total second virial coeffi- 
cient. At 500°K, it is 24 percent and at 1000°K 
it is 28 percent. At 400°K, the term x*Z12(E,/kT) 
which contains the 16th power of the dipole 
moment still contributes 3.1 percent to the 
second virial but at 575°K it is reduced to 
0.6 percent. Because of the importance of the 
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Fic. 2. Reduced second virial coefficients as a function 
of reduced temperature. The CC, CCC, --- are for the 
saturated hydrocarbons ethane, propane, etc.; C=C is 
ethene and CC=C is propene. 


dipole terms, the second virial coefficients really 
provide a sensitive method for the determination 
of dipole moments. 


II. GENERAL PROPERTIES OF SECOND VIRIAL 
COEFFICIENTS OF COMPLEX MOLECULES 


In making a general survey of the second 
virial coefficients of complex molecules, our first 
question is how well does the equation of 
corresponding states hold. In Fig. 2 we have 
plotted the reduced second virial coefficients, 
B(T)/V., of twelve substances as a function of 
their reduced temperatures, T7/T., where V, and 
T. are the critical volume and temperature, 
respectively. The data for some of the freons, 
as for example CCl:F2 (where the temperature 
of the measurements is uncertain to one degree) 
may not be very accurate. However, -a number 
of qualitative conclusions may be drawn: 

1. The equation of corresponding state holds 
well for temperatures above the critical. Below 
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the critical temperature there is a considerable 
spread in the values of the reduced second 
Virials. 

2. Molecules with large dipole moments have 
unusually large reduced second virial coefficients. 
However, from our previous section we see that 
the important quantity is not the dipole moment 
itself but x=y'E,,—*ro-®. Since E,, is proportional 
to T. and 7° is proportional to V., we should 
expect that the important reduced dipole 
variable should be 


2=(u'/T2V.7)10®. (8) 


Thus, the reduced second virial coefficients of 
H.0, NH3, and CH;Cl do not fall in the order 
of their dipole moments: 1.83X10-'® e.s.u., 
1.47X10-'8 es.u., and 1.89X10-'® es.u., re- 
spectively. But rather, they occur in the order 
of their values of z: z,0=1.32, 2yn,=.54, and 
Zcu,ci=.-40. 

3. The shape of the molecules does not greatly 
affect the reduced second virial coefficients until 
the temperature becomes considerably lower 
than the critical point. The curve marked 
‘Lennard-Jones spherical non-polar gas’’ is the 
reduced second virials for Ne, A, No, and CO 


BT) rennert-Jones = 0.68 Fi2(0.77T./ T). (9) 
spherical non-polar 

Deviations from this curve may be attributed 
to lack of symmetry of the molecules, the 
presence of dipole moments, or any other 
deviations from the form which Lennard-Jones 
assumed for the intermolecular forces. Actually 
the hydrocarbons and only slightly polar sub- 
stances have reduced second virial coefficients 
which lie more nearly along the dotted curve: 


B(T)/V.=0.89[ Fi2(0.73T ./T) 
— 0.00422 exp (4.47./T)]. (10) 


The development of the previous section makes 
it clear that the use of just the first power of 2 
in the second virial cannot hold at low tempera- 
tures for molecules with large dipoles. Actually, 
the deviation of individual molecules from these 
predictions usually exceeds the experimental 
errors so that it is necessary to fit each molecule 
as a special case. 

There are considerable experimental data for 
the olefin and paraffin series of hydrocarbons 
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which we can use to obtain additional generalities. 
For example, by comparing the virials for the 
four butene isomers which were studied by 
Roper, we see that the shape of the molecule 
does not greatly affect the second virials: 


B(225°K) B(280°K) B(335°K) 


cc/mole cc/mole cc/mole 
2-methyl propene —1436 —782 — 483 
1-butene — 1432 —777 —478 (11) 
trans-2-butene — 1504 — 854 —551 
cis-2-butene — 1645 — 868 —517 


Apparently the second virial coefficients for 
olefins are almost the same as for the corre- 
sponding paraffin. This is shown by a comparison 
at T=335°K: 

B(ethene) — B(ethane) = —118+149=31 cc/mole, 


B(propene)— B(propane) = —274+317=43, (12) 
B(1-butene) — B(n-butane) = —478+516=38. 


Perhaps the most curious hydrocarbon regularity 
is that the second virials are proportional to the 
square of the number of carbon atoms. Strict 
proportionality would be in contradiction to the 
theorem of corresponding states and seems quite 
unlikely. This is shown in Table III. All of these 
regularities seem to hold better at the higher 
temperatures. Before attempting an explanation, 
we should have experimental data for the higher 
hydrocarbons and for temperatures above the 
critical so that we can be more positive about 
our facts. 


III. THE IMPERFECT GAS AS AN 
ASSOCIATION PROBLEM 


An imperfect gas may be considered to be an 
equilibrium mixture of readily interconvertible 
single, double, triple, etc. molecules. Fortunately, 


TABLE III. 
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1. Paraffins 
n B(448°K) /n? B(523°K) /n? 
Ethane 2 —21.6 —15.0 
Propane 3 —17.8 —12.1 
n-Butane 4 —18.0 —12.4 
n-Heptane 7 —17.9* —12.9* 
2. Olefins 
n +B(225°K)/n? +B(280°K)/n? +B(335°K) /n2 
Ethene 2 —62.3 —41.6 —29.6 
Propene 3 —72.3 —44.1 —30.4 
2-Methyl- 

Propene 4 —89.7 —48.9 —30.2 
1-Butene 4 —89.5 —48.6 —29.9 
trans-2- 

— 4 —94.0 —53.4 —34.4 
cis-2- 

Butene 4 —102.8 —54.2 —32.3 








* Extrapolated from 548°K. 














Fic. 3. An imperfect gas as a mixture of single and double 
molecules (or collision complexes). 


only the double molecules can affect the second 
virial coefficient. In order to define the double 
molecules unambiguously, it is necessary to 
postulate that there exists a range of interactions 
such that whenever two molecules lie within this 
range, the energy of interaction is appreciable 
compared to kT and whenever they lie outside 
of this range their energy of interaction is 
negligible. If a molecule comes within the 
interaction range of another molecule the pair 
is considered to form a double molecule. This 
is illustrated in Fig. 3. The single molecules 
behave like molecules of a perfect gas until they 
collide with other molecules, whereupon we treat 
them separately as double molecules. The free 
volume in which the single molecules move like 
a perfect gas is V—b, where V is the molar 
volume of the system and )b, differs from the 
usual van der Waals 0 in having the range of 
interaction of the attractive forces substituted 
for the smaller collision diameter. Subtracting 
b, from V corresponds to the fact that single 
molecules are excluded from the definitional 
domain of the double molecules. It is true that 
double molecules behave like a perfect gas in a 
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volume slightly different from V—d, but this 
difference can affect only the third and higher 
virial coefficients. From these considerations, 
we see that up to the third virial coefficient, 
the equation of state of an imperfect gas may 
be written : 


p(V—ba) =(mi+n2)RT, (13) 


where 7; is the number of moles of single mole- 
cules and m2 is the number of moles of double 
molecules. The numbers of single and double 
molecules are determined by an ordinary equi- 
librium constant for association 


In the limit that V is very large, m; approaches 
unity and m2 approaches the value K(T)/V. 
Rearranging the equation of state (13) into the 
virial form, it is evident that 


B(T)=b.—K(T). (15) 


It is easy to derive this expression for the 
second virial coefficient directly from the 
fundamental statistical mechanical relation : 


N 
167? 





B(T)= f+ fa-exp (-B/er)yaaar. 


While this integration is to be carried out over 
the whole range of the variables, the predominant 
part of the integral comes from that range over 
which the molecules have an appreciable effect 
on one another. This is the same interaction 
range which we were discussing in the previous 
paragraph. Thus: 


N 
B(T) =—— [| --: | ((—exp(—E/kT))dQdr 
1) = f-- fa—exp (- Bye) 


interaction 
range 


=b,—K(T), (16) 


N 
b,= fo fasar (17) 
167? 


interaction 
range 


where now 





and 


K(T)= 





N 
ff exp (- Byer yaaa. (18) 
167? , 


interaction 
range 


It is evident that b, and K(T) have the physical 
significance which we attributed to them in the 


above discussion. The equilibrium constant 
K(T) for the formation of double molecules 
is the ratio of the partition function for 
double molecules to the square of the partition 
functions for the single molecules. Since very 
few of the double molecules exist in discrete 
vibrational or rotational levels, we may treat 
the problem classically. In the classical approxi- 
mation, it is evident that the internal degrees of 
freedom, the rotation of the double molecules, 
and the translational kinetic energy parts of the 
partition functions all cancel out of the equi- 
librium constant, leaving us with just the part 
of the phase integral for double molecules which 
is concerned with the translational potential 
energy. The additional factor of } arises from 
the symmetry number of the double molecule. 
Thus b, and K(T) are just the expressions to 
which our simplified discussion leads. 

Since K(7) is an ordinary chemical equilibrium 
constant, we can also express it in terms of the 
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Fic. 4. The “box” potential, an idealization of the 
energy of interaction of two molecules E as a function of 
their internuclear separation r. 


entropy AS and the energy AE of formation of 
the double from the single molecules. 


K(T)=exp (AS/R—AE/RT). (19) 


In general, AS and AE will depend somewhat on 
the temperature, but for comparatively small 
temperature ranges they may be considered to 
be constants. If we let b.c=exp (AS/R), the 
dependence of B(7) on temperature assumes the 
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functional form: 
B(T)=6,(1—c exp (—AE/RT)). (20) 


Here ba, c, and AE may be considered as three 
constants to be evaluated from the experi- 
mentally observed temperature variation of 
B(T). 

The physical significance of these constants is 
clarified by considering a hypothetical ‘“‘box”’ 
potential for which we can calculate B(T) in 
simple closed form. Here we assume that the 
molecules are spherically symmetrical; there is 
no energy of attraction until the molecules come 
a distance 7, apart; between 7, and the collision 
diameter 7» there is a constant energy of attrac- 
tion D; and it is impossible for the molecules to 
get any closer together than ro. Figure 4 makes 
this clear. Substituting this potential into 
Eq. (5), we get: 


2xNr,' %_°—To° 
B(T)= F [ -( —— } exp wer) (21) 
wT." 








Comparing this with Eq. (15) we see that: 
ba=2arNr,?/3, c=(ra?—ro*)/r, AE=—D. (22) 


Here we see once more that 0, is the van der 
Waals constant corresponding to the range of 
the attractive forces r. and AE is the energy of 
formation of the double molecules. However, 
a consideration of this box potential brings out 
the significance of c; it is the ratio of the free 
volume vy within the definitional domain of the 
double molecule to the total volume of the 
double molecule. Looking at these constants 
from a slightly different angle; bac=v; and so 
AS=R In vs, an expression similar to that which 
has become familiar in the theory of liquids. 
This interpretation of c in terms of a free volume 
is a general concept which should hold for any 
type of potential well. In case the molecule is 
asymmetrical, hindered rotation of the two 
molecules with respect to one another within 
the definitional domain of the double molecules 
will decrease the free volume and c. The extent 
of this decrease might be indicated by a “‘free- 
angle’ factor equal to the average solid angle 
through which the molecules can rotate in the 
double molecule divided by 4x. This would be 
analogous to the treatment of Kincaid and 


Eyring” of the theory of liquids made up of 
asymmetrical molecules. We must be cautious 
in our interpretations of c, as unusually small or 
unusually large values may be due to force fields 
of a peculiar nature as well as to shape factors. 

The exact second virial coefficient, B(T) 
=boF y(En/kT), for the spherical non-polar 
molecules which satisfy the Lennard-Jones 
potential (4) can be curve fitted to approximate 
B(T)=0b,(1—c exp (—AE/RT)). If we take bd, 
=4.75bo, c=0.835, and AE/R=—0.608E,,/k, 
the fit is perfect for E,,/kT =0.4988, 0.7046, and 
0.9953 and is better than one part in one hundred 
over the whole useful temperature range. Going 
back to the equation of corresponding states, 
we. could write these constants: 6,=3.23V.,, 
c=0.835, and E/R=0.463T./T. The constant 
c=0.835 is characteristic of the Lennard-Jones 
potential and complex molecules may have 
values quite different. Because of this variation 
in ¢ we must depart from the equation of corre- 
sponding states when we come to fitting the 
experimental values to within their experimental 
error. 

Unfortunately there are only a few molecules 
for which the experimental data are sufficiently 
accurate and cover a sufficiently large tempera- 
ture range so that three constants 0}, c, and 
AE/R can be uniquely determined. For example 
in the case of ethane, we might take: 


b,=135 cc/mole, c=0.731, 
and —AE/R=356 deg. 
or 


b,=470 cc/mole, c=0.841, 
and —AE/R=153 deg. 


Either set of constants would fit the experimental 
data to within five percent for a two-hundred- 
degree temperature range. It is therefore neces- 
sary to legislate arbitrarily the value of one of 
the constants, say b,. Making the most careful 
fits of the hydrocarbons and freons, we find that 
the average value of 6,=1.75V,. The value of 
ba=3.23V. of the Lennard-Jones gases is 
definitely too large and would lead to rather 
large deviations from the experimental values. 
This may not be apparent from our graph for 





J. Kincaid and H. Eyring, J. Chem. Phys. 6, 620 
(1938). 
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TABLE IV. 
Fraction 
double 
molecules 
at 1 
-—AE/R_ atmos. Ref. 
ba degrees and expt. 
cc/mole c K 298°K = data 
Ethane 251 0.778 244 0.018 6 
Propane 341 .681 347 .030 5 
n-Butane 452 .689 387 .047 7 
n-Heptane 749 559 629 141 8 
Ethene 223 .788 222 015 1 
ropene 316 .679 339 .027 1 
Propadiene 300* .614 382 .027 1 
1-Butene 430t 487 492 .045 1 
2-Methyl-Propene 430+ 493 490 045 1 
trans-2-Butene 430t .569 465 .047 1 
cis-2-Butene 430t 443 537 .047 1 
Freon—11 434 .729 399 .049 h 
CC1;F 
Freon—21 345 .920 306 .036 i 
CHCI1.F 
Freon—12 381 .631 345 031 j 
CC1.F 2 
Freon— 113 569 .888 335 .064 k 
CC1.F —CCIF2 
Methyl Chloride 239 582 469 .027 g 
CH;Cl 
Ammonia 62.8f .509 692 .013 f 
Steam 38.5 .420 1260 .045 3 








* Estimated from propane and propene. 

+ Estimated from ethene and propene. 

tba was varied so as to fit the experimental data best, all other 
values of ba in the table take ba =1.75Vc. 


the reduced second virial coefficient because of 
the small scale in which the points are plotted. 
By assuming a definite value of b,, the other two 
constants may be uniquely determined and have 
significance. 

In Table IV we have tabulated the best 
values of the constants b,, c, and AE/R for the 
hydrocarbons, freons, ammonia, steam, and 
methyl chloride. Let us interpret these results: 


1. Paraffins 


(a). AE/R increases with increasing number 
of carbon atoms. This corresponds to the long 
range of the van der Waals forces which enable 
distant atoms in a molecule appreciably to 
contribute to the energy of attraction of a 
molecule during collision. 

(b). The value of c decreases with increasing 
number of carbon atoms. This agrees with our 
expectation that c should decrease as the 
molecule becomes more unsymmetrical. 

(c). At one atmosphere pressure and room 
temperature, the fraction of double molecules 
increases with the number of carbon atoms. 
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AND WEEKS 


This may be considered a tendency for the 
higher hydrocarbons to cluster in the vapor 
phase. However, the hydrocarbons very nearly 
satisfy the theorem of corresponding states, so 
that this tendency would be the same if we 
compared the hydrocarbons at the same value 
of their reduced temperature, 7/7. 


2. Olefins 


(a). The value of AE/R increases and c 
decreases as we increase the number of carbon 
atoms. This is similar to the paraffin series. 

(b). cis-2-butene has a small dipole moment 
and its value of c is less than for the others. 


3. Freons 


(a). All of the freons have large values of c 
and approximately the same values of AE/R. 

(b). Because of the comparatively large errors 
in the experimental values, it was impossible to 
choose the values of c and of AE/R such that all 
of the experimental points were fitted to within 
5 percent. For the other substances tabulated, 
the maximum errors in the curve fitting were of 
the order of 1 percent. 


4. The Dipole Molecules: Steam, Ammonia, 
and Methyl Chloride 


(a). Methyl chloride has a large value of 
AE/R for a small molecule but otherwise its 
constants appear normal. Because it has such a 
small reduced dipole moment (even though its 
actual dipole is larger than that of water) it 
shows no extreme dipole behavior. 

(b). The values of 6, for steam and _ for 
ammonia are less than their critical volumes or 
about one-half of the normal values. If we had 
used larger values of b., we would have obtained 
a considerable deviation from thc experimental 
observations. 

(c). The values of AE/R for steam and 
ammonia are very large for such small molecules. 

(d). The values of c for steam and ammonia 
are abnormally small. This is to be expected 
since small values of c are a manifestation of 
hindered rotation. 

The large number of instantaneous double 
molecules in an imperfect gas is an interesting 
feature of our viewpoint. The fact that 4.5 
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percent of butane or 14 percent of heptane 
molecules at room temperature and one atmos- 
phere pressure in the vapor phase are in a state 
of collision should have some interesting conse- 
quences. We should expect deviations from 
Beer’s law for certain frequencies. Also, as 
H. Beutler and A. Turkevitch have suggested to 
us in a private discussion, some of the spectro- 
scopic absorption bands which are present in the 
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liquid but lacking in the molecular spectra 
should appear in the imperfect gas with an 
intensity proportional to the number of double 
molecules or proportional to the square of the 
pressure. 

We would like to express our appreciation to 
the Wisconsin Alumni Research Foundation for 
the financial support of one of us (I.F.W.) 
during the course of this work. 
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Temperature and Acidity 


W. B. PIrETENPOL 
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(Received December 11, 1941) 


Co-precipitation of radium and barium salts which ac- 
companies the cooling of a solution results in the deposition 
of a percentage of radium largely in excess of that of barium. 
Under similar conditions of cooling the relative amounts 
of radium and barium deposited have been determined at 
different temperatures throughout the cooling process and 
curves have been plotted. Determinations of the above 
nature have been made with radium-barium solutions 
containing different concentrations of hydrobromic acid, 
and the variation in the form of the curves is shown. The 
first crystals formed in the solution contain a larger per- 
centage of radium than succeeding crystals. It is pointed 
out that the efficiency of the refining process may be in- 


INTRODUCTION 


N the commercial refining of radium, the 

radium salts are commonly purified by frac- 
tional crystallization of a radium-barium solu- 
tion. The chlorides and bromides are usually 
employed and the details of the method have 
been fully described.! In a series of crystalliza- 
tions it has been pointed out? that the crystal 
fractions are progressively enriched with radium, 
while the liquor fractions are progressively im- 
poverished. The concentration of the radium at 
any step in the positive direction is expressed by 
the relation 

C=AK", 


where A is the initial concentration of the radium 


. Bull. U.S. Bur. Mines, No. 104 (1915). 
* Clarence E. Scholl, J. Am. Chem. Soc. 42, 889 (1925). 





creased by separating liquor and crystals at a point above 
room temperature. In this way the factor of enrichment 
is increased, thereby hastening the process. Factor of 
enrichment curves for solutions of different acidity show 
that for cooling from boiling to room temperature an 
increase in acidity decreases the crystallization factor. 
However the inverse relationship is shown to be true when 
the solution is cooled through a small temperature interval. 
The rate at which the solution cools is a factor in the 
relative amounts of radium and barium deposited. With 
rapid cooling a smaller percentage of radium is precipitated 
than with slow cooling. 


salt, m the number of crystallizations, and K the 
factor of enrichment or crystallization factor. K 
is defined as the ratio of the concentration of 
radium crystals separated to the concentration 
of radium in the original crystals. 

Nierman*® working with mesothorium and 
radium-barium bromide solutions, found that 
the crystallization factor is independent of the 
concentration of hydrobromic acid. He also 
concluded that for varying concentration of meso- 
thorium up to 2 mg per gram of salt, the crystal- 
lization factor remains constant. Unfortunately 
he made no reference to what precautions were 
taken in keeping the temperature limits constant, 
neither did he state that the rates of cooling of 
his solutions were uniform. As will be shown 


3 John L. Nierman, J. Phys. Chem. 24, 192 (1920). 
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later these are two important conditions to be 
considered in the determination of the factor of 
enrichment. 

Scholl? carrying out a fractional crystallization 
process with the radium-barium bromides and 
chlorides found that the value of K depended 
upon the temperature of maximum heating and 
minimum cooling, the speed of cooling, and the 
acidity of the solution. Scholl cooled his solutions 
saturated at 95°C by placing the crystallizing 
dishes for a period of 20 hours in a room of 
constant temperature at 20°C. He concluded 
that the greater the acidity the greater was the 
quantity of crystals formed and the smaller the 
factor of enrichment. Furthermore the effect of 
low concentrations of radium salt was found to 
be very marked if occurring with high concen- 
trations of acid. With these conditions it was 
indicated that the enrichment of radium is 
reversed, more radium occurring in the mother 
liquid than in the crystals. Scholl’s conclusions 
may now be interpreted in terms of the results 
presented in this paper. 


EXPERIMENTAL PROCEDURE 


A stock solution of radium-barium bromide 
was used, containing approximately 20 mg of 
radium per ton. 450 cc of this solution was placed 
in a beaker and heated in a bath the temperature 
of which could be controlled. This solution was 
heated to 95°C and under constant stirring was 
allowed to evaporate until saturation was ob- 
tained. The technique used to determine satura- 
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tion was to allow a slight stream of air to strike 
the surface of the solution and observe the 
formation of crystals. After many trial experi- 
ments this was found to be a satisfactory method 
of indicating saturation. 

At the point when crystallization was initiated 
samples were drawn by means of a graduated 
pipette to determine the radium-barium ratio. 
Because of the experimental difficulty in drawing 
two samples simultaneously, separate identical 
runs were made to determine the radium and the 
barium contents. For the radium determination, 
when the activity was high, proper dilution was 
made in order to obtain the correct amount for 
the most accurate measurement. The sample was 
then boiled gently to remove occluded radon gas, 
sealed in a flask under partial vacuum, and 
allowed to grow for a definite period. The 
amount of radium per cc of solution drawn was 
then measured in the usual way by an emanation 
electroscope of the Lind type. For the deter- 
mination of the barium in solution the samples 
drawn were evaporated to dryness under gentle 
heat and allowed to stand in a dry atmosphere 
for about two days before weighing. 

After the initial samples were drawn the 
solution was allowed to cool at a controlled rate. 
During the cooling the solution was constantly 
stirred and succeeding samples were drawn at 


spaced temperature intervals. However, preced- 
Curve Acidity 0.00 % 

‘ "6.75% 

0.45 % 
3.17% 
6.75 % 


Percent in Solution 


° 
100 90 80 70 60 $0 40 30 20 10 ° 


Temperoture, C° 


Fic. 2. 
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ing the drawing of a sample the stirring was 
temporarily discontinued to avoid the removal 
of minute crystals in suspension. Since the rate 
of cooling is a factor in the co-precipitation 
process, the cooling rate was reproduced as far 
as was possible in succeeding runs. The degree 
to which this was accomplished is shown by the 
curves in Fig. 1 covering the most critical portion 
of the temperature range. 

The acidity of the solution was changed by the 
addition of a certain percentage of hybromatic 
acid. While this gave an approximate value, it 
was not considered reliable because of the 
evaporation of water in initially reducing the 
solution to the point where it was saturated with 
barium bromide. A sample was therefore drawn 
at the end of each run by means of which the 
maximum acidity was determined. 


RESULTS 


The results of the experimental determination 
of the co-precipitation of radium and barium 
bromides as related to temperature and acidity 
are shown by curves in Fig. 2. In this figure, the 
percentages of radium and of barium remaining 
in solution are plotted as a function of tempera- 
ture. From these results the amounts of radium 
and barium precipitated in the form of crystals 
may readily be obtained and the factors of en- 
richment determined. For the various values of 
acidity in percentage by weight these factors are 
plotted in Fig. 3 as a function of the temperature 
in degrees C to which the solution was cooled. 

The conclusions to be drawn from the author’s 
results are obvious. In the first place, it is ap- 
parent that the factor of enrichment, as deter- 
mined by the radium-barium ratio in the crystal 
fractions, is greatly increased when the cooling 
process is stopped and the liquor is separated 
from the crystals at a temperature far above 
room temperature. Such procedure results in a 
far greater concentration ratio and, in addition, 
introduces a time saving factor, since under 
ordinary conditions the time rate of cooling is 
more rapid at high temperatures than in the 
neighborhood of room temperatures. 

The rate of formation of barium crystals is, 
however, a factor which cannot be neglected in 
the co-precipitation of radium. As an indication 
of this, the following experiments were per- 
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formed. A radium-barium solution was heated 
and evaporated to saturation as in the previous 
experiments. Samples were drawn to determine 
the initial radium and barium content. One-half 
of this solution was cooled rapidly by passing 
through a water cooler and the other half was 
cooled slowly as represented by the time curve 
in Fig. 1. After cooling through the same tem- 


perature interval samples were taken to deter- 


317%) by 
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Factor of Inrichment 


| 
0 90 80 70 60 80 40 30 20 


Temperature of Cooling 


Fic. 3. 


mine the radium-barium ratio and the acidity 
of the solution. In these tests the acidity was 
1.14 percent by weight. 


Test No. 1. Initial radium content at saturation, 
426.8 X 10-"° curie per cc. 
Factor of enrichment 


2.52 
3.29 


Temperature interval 


96.5° to 22°C 
96.5° to 22°C 


Cooled rapidly 
Cooled slowly 


Test No. 2. Initial radium content at saturation, 
451.8 107! curie per cc. 


Factor of enrichment 


2.66 
3.31. 


Temperature interval 


100° to 25°C 
100° to 25°C 


Cooled rapidly 
Cooled slowly 


The results show that with slow cooling a 
larger percentage of radium bromide is deposited 
with the barium bromide crystals. In terms of 
co-precipitation as an adsorption process these 
results are significant. However, in the commer- 
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cial refining process the variation in the crystal- 
lization factor with time rate of cooling may not 
be of great importance. 

The conclusions of previous experimenters 
working in the fractional crystallization of 
radium and barium salts may now be inter- 
preted in the light of the author’s results. Nier- 
man’s conclusions that the factor of enrichment 
is independent of the concentration of hydro- 
bromic acid and that the separation can be con- 
ducted advantageously in neutral or slightly acid 
solutions of the bromides, may be explained by 
reference to the curves of Fig. 3. Without careful 
control of the temperature interval through 
which the solution has been cooled no fair con- 
clusion concerning the effect of acidity on the 
factor of enrichment may be drawn. Further- 
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more, Scholl’s statement that “the greater the 
acidity the greater the quantity of crystals and 
the smaller the factor’ is confirmed by the 
present results for the special case, as in his 
experiments, when the solution was cooled from 
95°C to room temperature at 20°C. It may be 
mentioned that in the author’s experiments the 
concentration ratio is constantly changing with 
change in temperature. Work on the co-precipi- 
tation of radium and barium salts as affected by 
temperature and concentration ratio, and the 
nature of the adsorption process is in progress. 

In conclusion the author wishes to express his 
deep appreciation to Dr. O. C. Lester who sug- 
gested the problem and gave very valuable 
assistance in the experimental technique in- 
volved. 
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Intramolecular Energy Transfer 


The Fluorescence of Complexes of Europium 


S. I. WEISSMAN* 
Department of Chemistry, University ef California, Berkeley, California 
(Received December 17, 1941) 


The characteristic line fluorescence of trivalent europium is excited in certain organo- 
europium compounds by irradiation with light absorbed only by the organic part of the com- 
pound. The efficiency of excitation varies greatly with the nature of the compound, tem- 
perature, and solvent. Under optimum conditions, i.e., solution of a covalent compound at 
liquid-air temperature, almost unit efficiency has been obtained. The decay time of the 
fluorescence is independent of quenching; quenching prevents excitation of the europium ion. 
There seems to be a steady gradation of efficiency of transfer from covalent to ionic compounds. 


BSORPTION of energy by the trivalent 

ions of certain of the rare earths is often 
followed by emission of characteristic line 
luminescence. In most instances of rare earth 
fluorescence hitherto reported, the excitation has 
been accomplished by illumination of the 
material with light containing frequencies which 
are absorbed directly by the rare earth ions.! 


* National Research Council Fellow, 1941-1942. 

10. Deutschbein and R. Tomascheck, Ann. d. Physik 29, 
311 (1937); H. Gobrecht, Ann. d. Physik 28, 673 (1937). 
Except for the continuous absorption at wave-lengths lower 
than 2500A, the rare earth ion absorptions are weak and 
narrow, and as a result excitation of moderate intensities of 
luminescence require great intensities of exciting light. 


As we shall see, the excitation may be accom- 
plished, under suitable conditions, through light 
absorption by other constituents of the rare 
earth compound with subsequent transfer of 
energy to the rare earth ion. Indeed, in some 
cases excitation with wave-lengths in the near 
ultraviolet and the visible violet yields ‘‘for- 
bidden”’ rare earth fluorescence as bright as that 
of fluorescein under similar excitation. 

Experiments with compounds of europium 
only will be discussed here, although observations 
with samarium and terbium have shown that 
the phenomena occur in compounds of these 
substances also. 
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Comparison of the fluorescence of crystals of 
europium salicylaldehyde? and of europium chlo- 
ride hexahydrate illustrate the phenomenon. A 
monochromatic (30-40A wide) beam of light 
provided by a continuous source (the sun) and a 
Hilger quartz monochromator were focused first 
on a mass of crystals of EuCl;-6H.O. As the 
dial of the monochromotor was turned through 
the 3000A-—4000A region, bursts of the character- 
istic red orange fluorescence of europium* 
appeared at wave-lengths of exciting light 3920A 
and 3200A. These are the locations of the only 
fairly strong absorptions of europium ion in 
this region. When the europium chloride was 
replaced by a mass of crystals of europium 
salicylaldehyde, a compound possessing an 
intense and continuous absorption band between 
3200A and 4400A, characteristic europium 
fluorescence of almost uniform intensity appeared 
at all wave-lengths of exciting light between 
these limits. Although sharp localization of the 
act responsible for this intense absorption may 
not be possible, the internal (4f)® electronic 
system of the europium is certainly not involved 
since the same band appears in the salicylalde- 
hydes of colorless ions such as lanthanum and 
gadolinium. Hence it seems certain that excita- 
tion of the internal electronic system of the 
europium is provoked by light absorption in a 
region external to it. 

The yield of fluorescence under this mode of 
excitation, in all cases thus far studied, has been 
strongly dependent on temperature and, in 
solution, on the nature of the solvent. The 
intensity of fluorescence from a 5X10- molar 
benzene solution of europium salicylaldehyde 
increases by a factor of 5.1 as the solution is 
cooled from 45°C to 6°C.4 A solution of europium 
nitrate in water shows a change of only a few 
percent in its fluorescence intensity under the 
Same temperature change. The fluorescence 
yield from solid europium salicylaldehyde in- 
creases about 150-fold as the material is cooled 





* Eu(CsH,CHOO);, a non-ionic compound, insoluble in 
water, very soluble in benzene. 
*A line at 5790A, a group at 5880A, a group at 6100A. 
The measurements were made with a Weston photronic 
cell Model 594, type 2. The solution was contained in a 
yrex cell 1 cm deep; in such a cell there is total absorption 
of the exciting radiation—the 3200A-4000A region from a 
igh pressure mercury arc. 





OF COMPLEXES OF 





EUROPIUM 215 
from 300°K to 90°K; the yield from crystalline 
europium chloride hexahydrate increases by a 
factor of 2 at most over the same temperature 
range. 

The sensitivity of the fluorescence yield to the 
type of compound in which the transfer occurs 
may be seen from the data of Table I; here are 
given the intensities (in arbitrary units) at 
different temperatures of the fluorescence from 
2X10- molar solutions in benzene of europium 


TABLE I. Temperature dependence of fluorescence. 


1°¢ qi Ie I/I 
6.5 233 46 5.1 
19.0 183 26 7.0 
32.0 106 14 7.6 


i 


1, =Galvanometer deflection from fluorescence of europium nitro- 
benzoylacetonate. ; . 
I2=Deflection from europium salicylaldehyde. 


metanitrobenzoylacetonate and europium sal- 
cylaldehyde. The measurements were made 
under comparable conditions of total absorption 
of the 3650A radiation of a mercury arc. It will 
be noted that in. the temperature range studied 
the fluorescence from the nitrobenzoylacetonate 
is brighter than that from the salicylaldehyde 
but that the salicylaldehyde fluorescence in- 
creases more rapidly with decreasing temperature. 

Further, the yields from solutions of a given 
compound vary greatly from solvent to solvent. 
The fluorescence at room temperature from a 
10-* molar alcohol solution of europium benzoyl- 
acetonate, excited by radiation from a high 
pressure mercury arc, is barely visible in the 
dark, while the same material in benzene or 
toluene exhibits a fluorescence easily visible in a 
fully lighted room. However, the fluorescence 
from the alcohol solution increases very rapidly 
with decreasing temperature and at —80°C 
becomes as bright as the fluorescence from the 
toluene solution. Indeed, the fluorescence in- 
tensities from solutions at low temperatures of 
many compounds seem roughly equal. We would 
naturally suspect that in these cases unit quan- 
tum efficiency is approached. Measurement of 
the yield of the brightest fluorescence yet 
encountered supports this suspicion. The meas- 
urement was made by comparison with the 
intensity of fluorescence from a_ fluorescein 
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TABLE II. Relative efficiencies of fluorescence of various compounds.t 








Compound 


Formula of negative constituent 


Soln. in Soln. in 
a. e. i.* a. e. i.* 
room liquid 

temp. air 


Benzene 
solution 
room 
temp. 


Solid Solid 
room liquid 
temp. air 





Salicylaldehyde 


Benzoylacetonate 
Metanitrobenzoylacetonate 
Dibenzoylmethide 


—O—¢—COCH:; 


o-hydroxy acetophenone 


5 nitrosalicyl aldehyde 


3 nitrosalicyl aldehyde 


3 nitrosalicylate 


5 nitrosalicylate 


(NO2)s¢—O— 
NO2g—O— 
NO:g—O— 
NO»zg—O— 
(NO2)3¢—CO.— 
Pt(CN)6~ 


picrate 
o-nitrophenolate 
m-nitrophenolate 
p-nitrophenolate 

2-4-6 trinitro benzoate 


Platinocyanide 





¢—CO-=CH—CO—CH; 
NOsg—CO-=CH—CO—CH; 
¢—CO-=CH—CO—¢ 





a =not isolated as a solid. ‘ 
* The alcohol-ether-isopentane mixture. 


+ Details concerning the spectra and preparation of the individual compounds will be given in another communication. 


solution.® A 2.8X10-‘ molar solution of europium 
metanitrobenzoylacetonate in a 2-5-5 mixture 
by volume of ethyl alcohol, diethyl ether, and 
isopentane’ illuminated at liquid-air temperature 


5S. I. Vavilov has made measurements of the absolute 
vield from fluorescein solutions. S. I. Vavilov, Zeits. f. 
Physik 22, 266 (1924). 

6 This solvent has been extensively used in this labo- 
ratory by Professor Lewis and Dr. Lipkin. See G. N. Lewis, 
5 — and D. Lipkin, J. Am. Chem. Soc. 62, 2973 
(1940). 


by the 3650A line of a mercury arc yielded 
fluorescence 0.86 as intense as that given by a 
6X10-! molar solution of fluorescein in alcohol. 
Total absorption of the exciting radiation took 
place in each case. If the quantum yield of 
fluorescence of the fluorescein solution is 0.9, 
the yield for the europium solution is 0.85.’ 


7This is a quantum efficiency; account has been 
taken of the frequencies of the fluorescein and europium 
luminescences. 
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In those situations in which the quantum 
yield is less than unity, it is desirable to locate 
the fluorescence quenching act in the sequence 
of events after the initial absorption of light. 
Measurements of the variation of the lifetime of 
fluorescence with efficiency provide information 
on this point. The transitions responsible for the 
europium fluorescence are of the forbidden type, 
and we would expect an afterglow of duration 
10-*-10~? sec. Vavilov and Sevchenko* have 
found for aqueous solutions of inorganic euro- 
pium salts 7 (the decay constant) =0.84X10-* 
second. With a Becquerel phosphoroscope, I 
have found for the organic europium compounds 
described here values of 7 from 0.7X10-* to 
1.0X10- second. 

The phosphoroscope consisted of two 25-cm 
disks each containing four symmetrically spaced 
1-cm wide slots. The disks were driven at thirty 
revolutions per second by a synchronous motor 
which maintained them in exact synchronism 
with the flashes from a high pressure mercury 
arc. This synchronism is necessary so that the 
illuminating slots encounter the light beam at 
the same position in each cycle. For detection 
of changes in decay time—the point in which 
we are especially interested—two settings of the 
disks were used; at one setting the slots were 
staggered by 4.2X10-‘ sec., at the other by 
18 X 10-4 sec. The relative intensities of afterglow 
at these times were measured with a photronic 
cell and galvanometer. 

Following are a set of intensity readings taken 
on 1.4X10~ molar solution of europium meta- 
nitrobenzoylacetonate in the ether-alcohol-iso- 
pentane mixture at —50°C and at —190°C. 

t=4.2X10~ sec. t=18X10~ sec. 


— 50°C 8.8 , 2.2 
—190°C 19.5 5.0 


The ratios of the intensities at the beginning 
and end of an interval 1.4X10-* sec. long (two 
half-lives) differ by only two percent, while 
there is over a twofold change in efficiency at 
the two temperatures. 


*S. I. Vavilov and A. N. Sevchenko, Comptes rendus 
acad. sci. U. R. S. S. 27, 541 (1940). 


Similar experiments on europium picrate 
solutions in the alcohol-ether-isopentane solvent 
in which the intensities were compared by a 
visual polarization photometer, showed that as 
the efficiency changed by over fifty-fold the 
decay constant changed no more than thirty 
percent. These results indicate that the major 
part at least of the quenching of fluorescence is 
accomplished by some dissipative process oc- 
curring previous to the excitation of the rare 
earth ion. There seems to be remarkably little 
deactivation by non-radiation processes of those 
rare earth ions which do receive excitation 
energy. 

The dependence of the efficiency of energy 
transfer on the type of compound in which the 
europium ion and absorbing group are linked 
has already been mentioned. I have examined 
the fluorescence of the europium compounds 
listed in Table II. Estimates are given of the 
efficiencies of transfer of the solids and approxi- 
mately 10-4 molar solutions in the alcohol-ether- 
isopentane mixture at room temperature and 
liquid-air temperature; estimates on benzene 
solutions, where it has been possible to prepare 
them, are also included. The entire zero to 
unity range of efficiencies has been arbitrarily 
divided into four regions; fluorescence yields of 
twenty-five percent or higher are labeled 3, zero 
yield 0, and intermediate yields 2 and 1. 

A steady decrease in efficiency of energy 
transfer from the covalent, benzene soluble, 
benzoylacetonates, to the ionic crystals of the 
platinocyanide is apparent. Evidently energy 
does not easily cross ionic bonds. Just what 
conditions of conjugation and resonance are 
necessary for ready transfer across non-ionic 
bonds has not yet been established. Further 
work on this point is in progress. 


ACKNOWLEDGMENTS 


All the europium used in this research was 
supplied by Dr. H. N. McCoy. A portion had 
been originally given to Dr. Simon Freed from 
whom I borrowed it, and the rest was given me 
by Dr. McCoy. It is a pleasure to acknowledge 
the generosity of Drs. McCoy and Freed. 











APRIL, 1942 


JOURNAL OF CHEMICAL PHYSICS 


Frequency Spectrum of Crystalline Solids 


ELLiott W. MontrROLL* 








VOLUME 10 


Sterling Chemistry Laboratory, Yale University,t{ New Haven, Connecticut 


' (Received January 2, 1942) 


It is generally accepted that a normal crystalline solid 
can be pictured at absolute zero as an assembly of mole- 
cules arranged at periodically placed lattice points. Since 
at higher temperatures each molecule becomes a harmonic 
oscillator about its lattice point, in order to calculate 
thermodynamic properties of the crystal it is necessary 
to know the distribution of its internal normal modes of 
vibration. On the basis of the Born-von Karman model 
these normal modes of vibration are roots of a secular 
determinant. In this paper it is shown that the 2th 
moment of the distribution function of normal modes is 
proportional to the trace of the mth power of the matrix of 
the Born-von Karman determinant. By expressing the dis- 


I. INTRODUCTION 


O* of the failures of classical mechanics and 
classical statistics was their inability to 
explain even qualitatively the decrease in specific 
heat of solids with decreasing temperatures; and 
one of the selling points of Planck’s original 
quantum theory had its source in the success of 
Einstein’s! application of Planck’s ‘‘quantized 
oscillators” to the lattice vibrations of a solid. 
By Einstein’s method the vibrational factor of 
the partition function of a normal crystal is 


3N exp (—hyv;/2kT) 
- \ 


i=1 1—exp (—hy;/kT) 





and the Einstein assumption of equal numerical 
value of all 3N frequencies of a crystal explains 
the observed decrease of the specific heat below 
the Dulong-Petit value of 3NkT at low tem- 
peratures. Here h=Planck’s constant, k= Boltz- 
mann’s constant, 7 = absolute temperature, 
v:=frequency of the ith normal mode of vibra- 
tion, N=number of lattice points. 

Refinements in specific heat measurements in 
the first decade of this century demonstrated 
that the agreement between experiments and 
Einstein’s theory was not quantitative at low 


* Sterling Research Fellow in Chemistry. 
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911) Einstein, Ann. d. Physik 22, 180 (1906); 34, 170 
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tribution function as a linear combination of Legendre 
polynomials it is shown that the coefficient of each poly- 
nomial is a linear combination of the moments. The fre- 
quency distribution function of a two-dimensional simple 
cubic lattice is calculated by the above method and turns 
out to have two maxima. Usually the equation for a ther- 
modynamic function F(T) involves the integral of the 
product of the frequency distribution function g(v) and 
a known function K(T, v). We show here that when F(7) 
is a known function of T an integral equation results with 
g(v) under the integral sign. This integral equation can be 
solved for g(v) by use of Fourier transforms. 





temperatures; and Debye? and Born and von 
Karman* pointed out that the discrepancies 
were probably due to the existence of a distribu- 
tion of normal modes of vibration instead of a 
single ‘Einstein frequency.’ Thus, for example, 
the temperature dependent vibrational contri- 
bution per lattice point to the free energy should 
be written! 


Fy(T)= -er f g(v) log (1—e-"”/*T)dp, (2) 
0 


where g(v)dv is the fraction of normal modes 
with frequencies between v and v+dy». 

Debye’s proposition to treat the solid as a 
continuum with the same elastic properties led 
to a frequency distribution proportional to v? and 
to equations of physical quantities in terms of a 
parameter 0 =hv,/k which theoretically should 
be a constant for a given crystal. The term v, is 
the largest frequency of the normal modes of 
vibration, and @ is usually referred to as the 
Debye or characteristic temperature. Unfor- 
tunately this development which mathematically 
is quite simple is insufficient, for experiments 


2 P. Debye, Ann. d. Physik 39, 789 (1912). P 
3M. Born and T. von Karman, Physik. Zeits. 13, 297 
(1912); 14, 15 (1913). 

‘For a general discussion of theory of normal solids cf. 
E. Schroedinger, Handbuch der Physik (1926), Vol. 10; 
M. Born and M. Géppert-Mayer, Handbuch der Physik 
(1933), Vol. 24; R. Fowler, Statistical Mechanics (1936), 
Chap. 4. 
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show that © is not constant over a large tem- 
perature range. 

The extension of Born and von Karman. of 
Einstein’s work afforded a more accurate 
physical picture than that suggested by Debye. 
However, it led to a mathematical problem 
which could be solved analytically for a one- 
dimensional solid but which required tedious 
numerical calculations in two- and three-dimen- 
sional crystals. In the Born-von Karman theory 
one sets up the equations of motion for a set of 
coupled oscillators with the same coupling and 
force constants as those of the crystal, and shows 
that the normal modes of vibration are the 
characteristic roots of a secular determinant of 
order N*. Actually the determinant can be 
factored into secular determinants of very small 
order, but to find g(v) by present methods the 
roots of all the resulting equations (and there 
are O(N*) of these) must be calculated and the 
number of roots in every small frequency interval 
must be determined. Even though it seems 
reasonable to base the g(v) computations on the 
solution of a large random sample of cubic 
equations, the lack of popularity of the Born-von 
Karman approach is understandable. 

The first part of this paper will be devoted to 
developing an analytical method of finding g(v) 
in a Born-von Karman crystal. This method will 
depend on the matrix theorem that the sum of 
the kth powers of the characteristic values of a 
matrix is equal to the trace of the kth power of 
the matrix. Noticing that the sum of the kth 
powers of the characteristic values of the Born- 
von Karman matrix is proportional to the kth 
moment of g(v), and that a continuous distribu- 
tion function defined in a finite interval can be 
expanded as a linear combination of Legendre 
polynomials, it is apparent that the coefficients 
of the Legendre polynomials can be expressed as 
linear combinations of the moments of g(v). This 
means that g(v) can be expressed analytically as 
a function of the lattice constants of the crystal. 

These observations will be applied in detail 
to a two-dimensional cubic lattice, postponing to 
some later date the application (which is not much 
more difficult) to real three-dimensional crystals. 

The second part of this paper will concern the 
determination of g(v) from experimental specific 
heat data by use of Fourier transforms. 





FREQUENCY SPECTRUM OF CRYSTALS 
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II. g(v) FOR THE BORN-VON KARMAN MODEL 
1. Mathematical Techniques 


The possible frequencies of normal modes of 
vibration in a Born-von Karman crystal will be 
shown in the next section to be roots of a secular 
determinant (see Eq. (17a)) and two approaches 
to the problem of deriving their distribution from 
the determinant are: (a) by calculating all or a 
large random sample of its roots and counting 
the number in each small frequency interval; 
(b) by evaluating the moments of the distribu- 
tion function and from the moments determining 
the distribution. 

Blackman’s®’ researches have demonstrated 
that the first method involves the solution of a 
large collection of quadratic or cubic equations 
and thus an enormous amount of tedious nu- 
merical work—and further, they show that the 
computation must be repeated for every set of 
lattice parameters, so method (a) cannot ezsily 
show the influence of changes in lattice param- 
eters on changes in physical properties. Here we 
shall devise a technique based on the second 
approach which avoids tedious numerical work. 

Let M be the matrix of a symmetrical secular 
determinant whose roots are {\;}. Then 


M¢ =) itt, (1) 


where ¢; is the ith characteristic vector of the 
matrix, and X; is the corresponding characteristic 
value. Multiplication of both sides of (1) by M 
yields 
Mt; = MA ab; =\;Mu; 

=)ith; 
and in general 

M*t;=),*t;. (2) 


Thus \;* is the characteristic value corresponding 
to the characteristic vector ; of M*. Now it is 
well known* that the trace of a matrix is the sum 
of its characteristic values, so 


Trace M‘=)}°; A+. (3) 
Since the total number of characteristic values 


(where a characteristic value of degree of de- 
generacy d is counted d times) is equal to the 





5M. Blackman, Proc. Roy. Soc. 148, 365, 384 (1935); 
159, 416 (1937) ; see also P. Fine, Phys. Rev. 56, 355 (1939). 

6 Cf. Rojansky, Introductory Quantum Mechanics (New 
York, 1938), p. 324. 
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order of the matrix, which we shall assume to be 
m, the average values of powers of the charac- 
teristic values of M are 


> A\*/m= (Trace M*)/m. (4) 


But this is exactly the kth moment of the dis- 
tribution function of the characteristic values 
and will be abbreviated by px. 

We must now find a way of expressing the dis- 
tribution function g(A) in terms of the set of 
moments {ux}. Let As and Az be, respectively, 
the smallest and largest characteristic values of 
M (since the matrix is finite they actually exist; 
also, since the matrix was assumed to be sym- 
metrical, the \’s are all real) and let us assume 
that the elements of the set {d;} are densely 
distributed. 

Suppose the distribution function is expanded 
as a linear combination of Legendre polynomials 


x 2A—(Azr+As) . 
gA)=d anPa( ), (5) 
n=0 AL—As 
1 @ 


P,(x) =——_ 
2"n! dx" 





where 





(a°—1)* (6a) 
and 


f Pr(x)Pj(xe)dx=26,;/(2k+1). (6b) 


To evaluate the coefficients {a,} multiply both 
sides of (5) by 


2A\—AL—As 
p,(——~—*)an 
Ar—As 


and integrate from As to Az 
AL 2\—AL—As 

f gayPi(——*—*)an 
As AL—As 


rd AL 2A\—AL—As 2A\—AL—As 
“inf CECA) 
n=0 AS AL—As AL—As 





Letting x=(2A—AL—As)/(AL—As) 
(Ar—As) 
“en g(A)Pi.(x)dx 
wail 





Ps 1(AL—A; 
=La f | - ») D(a) Pu(x)de 


=ax(Ar—As)/(2k+1). 
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Therefore 
(2k+1) f) ((it+x)ArL+(1—x)As 
a=—— f ef ) Paced: 
2 4 2 (7) 


If we define 


1 
Uy. -{ x*g(A)dx 
1 


2 AL 
nen f (2—A,—As)*g(A)dd 
(AL—As)*t! Sas 
2 AL 
-—— [ (polynomial in \)g(A)dX, (8) 
(A, —As)**! AS 


it is apparent that the a,’s can be written as 
linear combinations of the u,.’s. Furthermore, the 
definition of g(\) implies that its moments are 


AL 
‘ign f g(A)Nidd 


AS 
(9) 
= (Trace M‘)/m, 


so the u,’s are linear combinations of y;’s and 
therefore the coefficients a, are also linear com- 
binations of the moments (or the traces of M’‘) 
ui of g(Q). 

In the theory of frequency spectrum of solids 
the frequencies of normal modes are proportional 
to the square root of the roots of the secular 
equation, so the smallest frequency is the nega- 
tive of the largest, and the odd moments of the 
distribution function are zero. Considering this 
case in detail we have: 





g(A) =D anP(d/Ax) (10a) 
n=0 
with 
(2n+1) ¢' 
a,=— ff eer Patards. (10b) 
<i 
If 
1 
m= f gGr)atdx, (11) 
ai 
then 
AL 
m= f g(A)A‘dd 
—AL 
“=z ttly,, (12) 
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Therefore 


1 1 
oe=- f g(xX 7) ‘Idx 
2/7_; 


= po/ 2X1 


3 7 
a,=- f g(x) -xdx 
2/7_; 


3 
=~(ui/d12) =0. 
2 


Similarly, since the odd Legendre polynomials 
involve only odd powers of x, and since the odd 
moments of g(A) are zero: 


d)=A3="* ++ =d2n41=0, 


5 f' 3 1 
a2, =- f g(x.) (xt , a 
2/7_, 2 2 


= 5(3p2 ans MoA 1”) /4) “*, 
= 9( 354 — 302d 1? + 3y0dz') /16X25, 


A¢=13(231pr5— 315 psr22>+ 105 por! — 51°) /32A77, 








Miiy, mt a(t, m — U1+1, m+ U1, m— Uji, m) FY (U1, m +1, m 


— U 1-1, ae eet m—-1 + U1, a~tia U1+1, m+1 0141, m—1 FU, m—V1,m— uU 1-1, m+1 01-1, ee) _ 0 


and 


OF CRYSTALS 


1 d= 
f g(adn| — (x?—1 ya 
a’ dx** 


4k+1 — ; .| 
a ) 24 \ 
(2k) !27*+1 7 ldy** 


where after the differentiation is performed in 
the last expression, y* is replaced by u,/Xx*. 


and in general 


(13) 


k P . 
I, =ux/r7", 


2. Two-Dimensional Cubic Lattice 


Here we shall apply the Moment-Trace 
method to a_ two-dimensional one-component 
cubic lattice of the Born-von Karman type which 
includes interactions between nearest and next 
nearest neighbors. Let a be the binding force 
between two nearest neighbors whose distance 
of separation is a and let y be that for two next 
nearest neighbors whose distance of separation 
is av2. If there are 2N X2N particles of mass M 
in our square lattice, each atom can be specified 
by two letters /,m (1=l,m=2N) and its dis- 
placement components along rows and columns 
can be represented by 2; » and 2%, », respectively. 
The motion of the particle (/, m) is described by 
the simultaneous differential equations® 


— U1+1, m+1— V141, m+1 + ui, m+, m 


(14a) 


Mi, mta(v;, m— Vi+1,m +7, m— Vi-1, a) +y(u l,m +), ie U1+1, m+1— V1+1, m4+1 + U1, m+, m 


~*~ 35.2 ang Vi—-1, m—1 + U1,m— V1, m — U141, m—-1 +0141, m—1 + ii, m— V1,m— Vi-1, m+1 +v)-1, m+1) = 0. 


By assuming the periodic solutions 


Uim=u' exp i(2xvtt+1l¢itme.), 


Vi,m=v exp i(2rvt+lei+me¢s), 


¢gi=7a1/N; 


where a; and a2 are integers satisfying 


¢2=1d2/N, 


(14b) 


(15a) 
(15b) 
(15c) 


—N=a,, aN, 


and substituting (15) into (14) one can easily show that the frequencies of the normal modes of 
vibration of the crystal must be roots of the characteristic equations 


| A(¢1, ¢2)—40*Mv* 
| B(¢a, $2) 


B(¢2, ¢1) 


| =0, (16) 


A(¢2, ¢1)—4a°Mv? 


A(¢1, y2) =2a(1—cos g1) +4y(1—cos ¢1 cos ¢2), 
B(¢1, g2) =4y sin ¢1 sin ¢o. 
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Henceforth it will be assumed that v is measured in units of 27\/M so that the factor 422M can be 
omitted. If the matrix of the determinant (16) is represented by M?(a,, a2), the matrix of the entire 


characteristic determinant is 
0 


0 
L 0 





(M?(—N, —N) 


0 
M*(—N, —N+1) 
0 
0 


M?(N, N-1) 
0 


0 





M?(N, N)) 


(17a) 


(the submatrices are represented as squares because the roots of their characteristic determinants 
are proportional to the squares of the frequencies of the normal modes of vibration). Of course: 


(M**(—N, —N) 
0 


0 


0 
M*(—N, —N-+1) 


0 


0 
0 


M*(N, N—1) 





0 


\ 


0 0 


0 
0 


0 


\ 





M*(N, N)) 


Since all the small °(a;, a2) matrices are of the same form, namely 


M , a , a mat 
1 2 


A* 
(where A =A(¢i, ¢2); A*=A(¢2, ¢1) ; B=B(¢1, 2) =B(¢o, ¢1)), 
N N 
TraceM*= } > Trace M**(aj, ae) 
ai=—N ag=—N 
and we first need only to find traces for one general pair (a1, a2). Now 


Trace M?(a1, a2) =A(¢1, ¢2) +A(¢2, ¢1) 
Trace M?= >> [A(mai/N, ma2/N)+A(2a2/N, rai/N) } 


a1, a2 


and 


=2 >> A(xa,/N, ra2/N). 


a1, a2 


M*(a, @2) = M?(ay, a2) + M*(a4, a2) 


+ 
+ 


ve 


A?+B? 
B(A+A*) 


A B 

B ) 
B(A+4A*) 
B?+(A*)? 


) 


Trace M‘= >> [A?+B?+B?+(A*)?] 


a1, a2 


=2 > (A*+B?’). 


ai, a2 
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Continuation of this process leads to 


Trace M*=2 > (A*+3AB?), 


ai, a2 


Trace M§=2 >> (A‘+4A?B?+2B?AA*+B‘), 


a1, a2 


Trace M'°=2 > (A'+5A*B*+5A°B?A*+5A B‘). (20e) 


ai, a2 


These Traces will be sufficient to yield a g(v) correct to the 10th Legendre polynomial. In general 
the sum of the diagonal elements of the matrix 


A B A B A B 
M* (ay, a:)=( )( ): . ( ) 
B A*/\B_ A* B A* 


xX - 





k matrices 
is Trace M**(a, a2). 
Actually 
Trace M?=2 >> {2a(1—cos ¢;)+4y(1—cos ¢; cos ¢2) } 


a1, a2 
2{ ¥ (2a+4y)—4a © cos gi—8y(d cos ¢1)?}. 
a1, a2 ai, a2 a1 
But as shown in the appendix 
N 


> cos gi=—l. 
aij=—N 


Therefore 
Trace M? 2a 4y 


7 = (2a+4y)+ - 
2(2N+1)? (2N+1) (2N+1)? 





Since N is very large, and since this is exactly the second moment of g(v): 


Trace M? 
be =———— = (2a+ 47). 
2(2N+1)? 
To find ys we must first calculate 


Trace M*=2(>, A?+)>. B?). 


For this calculation we find 


> A?= ¥ [4a*(1—2 cos gi+cos* g2)+16ay(1—cos ¢1)(1—cos ¢1 cos ¢2) 
aj, a2 
+16y?(1—2 cos ¢; cos ¢g2+cos* ¢; COs* ¢2) }. 
But from the appendix: 


> cos? g=N+1, > cos ¢gi=—1, 


a1 


so 


> A*=4a2(2N+1)?+4a2(2N+1)(N+1)+16ay(2N+1)?+ 16y2(2N+1)?+ 16y2(N+1)?+0(N) 


@1, a2 
and 


Al LD A*/(2N+1)?=4(a+2y)?+2(a?+ 27’) +O(1/N). 
so, 


LD B?=16y7°(d sin? ¢1)’, 
a1 


a1, a2 
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but from the appendix a; sin? g;= N, therefore 
> B*/(2N+1)?=4y’. 





mer ps= (Trace M*)/2(2N+1)?,  ws=4(a+2y)?+2(a?+4y?)+0(1/N). (21b) 
In a similar manner some of the higher moments can be shown to be (to within terms of order 1/N) 
ue = 8(a+2y)*+12(a+2y)(a*?+4y’), (21c) 
us = 16(a+2y)'+48(a+ 2y)*(a?+4y") +6(a?+4y*)?+400*y’, (21d) 


Lio = 32(a+2y)>+ 160(a+ 2y)3(a?+4y") + 12077(a+2y)(7a?+8y") 


+40a*y?(a+y7)+60a'(a+2y). (21e) 
Using the results of the previous section, if we assume 


Dn 


g(v) = Lo denPon(v/ vz) (22) 


n=0 
(v_=largest frequency of a normal mode) the first few coefficients are 
ay=0.5/vz, 
a2=5(6a+12y— v1”) /40,3, 
a4=9[140(a+2y)?+70(a?+4y*) —60(a+2y)v1?+3v7' ]/16r7>. 
The number of parameters involved in the a’s can immediately be reduced by one, for as Blackman’ 


has shown, the largest normal mode v, results from the 2X2 determinant in which g:=0 and g2.=7 
(or vice versa, g2=7 and g2.=0); therefore by (16). 


vp>=4(at+2y). 
This means 
VviLag= 3, 


via2= 3 
vds=9(9a?—104ay+36y")/8(a+2y)?, etc. 


Since the purpose of this paper is more to introduce new methods than to obtain particular 
results, iet us limit ourselves to the special case y = 0.05a. This a, y relationship is of special interest 
because it is the one used by Blackman in his detailed calculations—and may therefore give us an 
opportunity to compare the results of the trace-moment approach with those obtained by Blackman 
in his detailed root calculations. In this case v,?=4.4a. 

The first six even moments in the special case y=0.05a are’ (of course the odd moments are all 


zero) 
bo= 1, pe = 0.28150826r,5, 
we = 0.577", ps =0.23560495 75, 
us =0.35433884y,', io =0.20316225r,"° 
and the corresponding Legendre coefficients of (22) are ) 


vpay=0.5, viag= 0.3704, 
v,d2=0.625, v1ag= —0.3182, 
v,a3=0.2260, Viai9= — 0.0133. l 


fay large number of significant figures is necessary because one must take differences of multiples of the y's in order 
to find the a’s. 
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The resulting function v ,g(v) is plotted in Fig. | | ] 
1. It qualitatively resembles the equivalent curve 
of Blackman as obtained by direct numerical com- 1.2,- 
putation of a random sample of normal modes. 
Blackman’s first peak is sharper and is displaced 
a little more to the left. The quantitative disagree- 


ment may be the result of Blackman’s smoothing a 7 
process, or it may have its source in the existence YOY) 

of a large value for one or more of the coefficients 

of the higher Legendre polynomials that have been 6 f- ~ 


neglected here. Since our g(v) is not quite zero at 

v=0 as it should be, we must be omitting some 

contribution of the higher polynomials. 3 
The calculation of physical quantities from a 

frequency distribution of the form (22) offers no 














special difficulty. For example the specific heat | | | 
is in general 0 25 5 15 10 
vl h*v*g(v)dv V/V, 
Cok f . . /¥, 
0 4k®T? sinh? hy/2kT Fic. 1. 


so in a low temperature region where one may expand (sinh? hv/2kT)-! as an exponential series, 
elementary integrals of the form 


VL 
f e—"’/kT (polynomial in v)dv 


0 


are all that are involved, while at higher temperatures where it would be desirable to use a power 
series expansion for (sinh* hvy/kT)—! one encounters only sums of integrals of the form 


VL 
f (polynomials in v)dv. 


0 


Before concluding this section it might be worth while reminding the reader that the expressing 
of g(v) as a linear combination of Legendre polynomials is somewhat arbitrary. Any other artifice 
which allows one to express g(v) as a function of the moments would be equally valid. The Legendre 
development merely seemed to be the most straightforward and simplest approach which occurred 
to the author. 


III. FREQUENCY DISTRIBUTION FROM EXPERIMENTAL DATA 


If in a crystal of N structureless atoms or in a complex molecule of N atoms, the frequency dis- 
tribution function g(v) is defined so that 3Ng(v)dv is the number of normal modes of vibration in the 
frequency interval vy and v+dy, the logarithm of the partition function (per normal mode) is given by 


log Z(T) = -{ g(v) log (1—e-’"/*") dv (23a) 


0 


under the restriction 


f g(v)dv=1. (23b) 
0 
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Formally a knowledge of Z(T) at all temperatures permits (23a) to be regarded as an integral 
equation in which g(v) occurs as the unknown function. Also, the average energy per normal mode is 


E(T) =kT?0 log Z/dT 


» hvg(v)d 
-{ ee. (24) 


elv/kT _ 1 





(using for zero energy the state of lowest energy of the system including the residual energy hv/2 
per normal mode) and by definition the specific heat is 


C,(T) =(0E/8T), 


hv \? g(v)dp 
-* ff ( —-) , (25) 
2kT/ sinh? (hv/2kT) 


Making the substitution @=h/kT and writing 
log Z(T)=2Z(0); E(@) 



















C(0) =C.(T)/k, 





=0E(T/h; 








we have the three fundamental integral equations 


a 


Z(0)= -{ g(v) log (1—e~*) dy, (26) 


0 






E()= { g(v)Ovdv/(e%”—1), (27) 
0 


* @y?g(v)dv 
C(e) = f 
> 4sinh2 (6/2) 











all of which are of the form 





400) = f g(v)K(@v)dv. (28) 


0 






It is quite clear that if the partition function, the internal energy or the specific heat of a given 
system could be found at all temperatures, and if the above integral equation could be inverted* 
so that g(v) would be a function of the measurable quantities it would be possible to compute g(v) 
and perhaps get a better understanding of the internal structure of a given crystal or complex 
molecule. 

As has been demonstrated in Appendix II, the formal solution of (28) is 








~ du 


eo)=— f f(0) (v0) "dé, (29) 
a T(u) 





where 


wef xK(x)dx. 


0 





Inasmuch as? 





BD 


f x™ log (1—e7*)dx = — ¢(2+7u)T'(1+7u) 


0 





* Inversions of the partition function have been made before by S. Bauer, J. Chem. Phys. 6, 403 (1938); 7, 1097 (1939). 
. These integrals are easily evaluated by expanding the transcendental functions in the integrands as sums of powers 
a ¢™. 














FREQUENCY SPECTRUM OF CRYSTALS 


(¢(z) is the Riemann zeta-function defined by ¢(z) =>" n-*) 


n=1 


© g@tdy »  gint2dy ; : 
f = ¢(2+iu)T(2+in), f ~= ¢(2+iu)T(3+7), 
oe &—! 0 ry sinh? x/2 


the formal equations for g(v) are 


1 Lt 
g(v)= [- ‘. Z(0)( v8)" 
2a » ¢(2+iu) P(A +in) 





1 * du 
g(v)= —{ fe 8) (v0)'"dé, 
2a J_, ((2+in)T (2+iu) 


1 e du 
g(v) = f f C(8)(v6)'"dé 
2nr J_, ¢(2+iu)P(3-+in) 


As is to be expected, substitution of the Debye specific heat formula 


12 L x3dx 36v, 
C(0)=C,(T) vK=3( J _ ) 
Ay, t,o et—1 eh a 


into (32) can be shown to yield the frequency distribution 


[992/01 if vine 


g(v) = (33) 


lo if »>»:, 


For actual calculation of g(v) from (32) there seem to be several possible approaches. One might 
integrate (32) numerically after substitution of the experimental C(@) values derived from specific 
heat measurements. However, since most specific heat curves are almost of the Debye form, one 
might divide C(@) into a Debye and a residual portion. The Debye portion would contribute the 
amount (33) to g(v) and the rest of g(v) might be found either by numerical integration or by ob- 
taining empirical equations for the residual C(@). In the latter case one would attempt to fit a power 
series in @ to the high temperature part and a power series in e~* in the low temperature region. Then 
32) could probably be integrated analytically. Another possibility would be to assume that g(v) 
could be expanded in Legendre polynomials and finding the coefficients by applying the orthogonality 
property of the polynomials to the assumed equation (32). 

In conclusion the author would like to thank Professor J. G. Kirkwood for his discussions con- 
cerning this problem. 





APPENDIX I. SUMMATION FORMULA 


[In the main part of this paper the following trigonometric To sum C; and 5S; it is convenient to find first 
sums are needed 
N N N 
C.= = cost =: S.= > ° E.= =. _exp (ikwa/N). 


a=-N re a=—! 4 a=— 








228 


If x=exp (ikr/N) #1 when 0<k<N, 


N 
Ey= 2D xt=x-N(1—x2N4H)/(1—x) 
a=—-N 
= e~itk( | — e2rikgike/N) /(1 — eiktiN), 
Since 
e~itk = (—1)4, 
E.= (- 1)* 


when O0<k<N. 


Obviously Ex5=2N+1. 
Now 


N 


O=% b (etta/N + e-ita/N)k 
a=—N 


k N 
1 k! 
on Sn 


eit (2m —k)/N, 
2*n=om'!(k — m) !a=—N 


When bk is odd, 2m—k is also odd and 


N 
= exp {in(2m—k)/N}=-—1. 


r= 


12 k! 


 aik—min 


Coad) re ae 


When & is even (2m—k) is also even and it is zero when 
k=2m; so 


we rm 4 k! 
m=0 mi(k—m)!" «in+im!(k—m)! 


Cceven) ad | 


4 (2N+1)k! 


(k/2)!(k/2)! 





ae ely 


~ 2k) —imi(k—m)! 


“ 


2Nk! 
(k/2)'(k/2)! 


Citecens = 1 +2!-*N(k!)/(k/2)!'(k/2)!. 


To find S; we proceed in essentially the same manner: 


N 
Sk= Oe D (eita/N —e-ina/N)k 


a=-N 


1 4 k! N seni nies 
“am, 2,alb—wyi,2 mr) 


- 1 > _ki(—1)" vy et ta(k—2m)/ N 
(21)*,-=9m!(k—m)!,—y , 


When bk is odd k—2m is odd and 
—1 * (—1)™k! 

Sr oac awn Ak c eae Spoememaee 
odd) (2i)*,~,m\(k—m)! 


=~ gpill- Dk=0. 
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MONTROLL 


If k is even k—2m is also even 


| 2, ki(—1)™ | & oRM—1)™ 


Sk ta Bo ne + ee, x | eee 7 ee fy 
feven) (21)*| ,.< m!(k—m)! "mace 41ml(k—m)! 





(—1)*/2(2N+1)kR! 
(k/2)'(k/2)! 


(2i)*| ~ym!(k—m)!" (k/2)'(k/2)! 


1 | & eu—1)m | a(2N)R! 
~ (2i)* 


7 7 N(k!) j 
~aveven) ™ Fk-1(k/2)1(k/2)! 





The first few C;.’s and S;’s are 


k Ck Si, 
0 2N+1 2N+ 
1 —1 0 
2 N+1 N 
3, 5,4, 9 —1 0 
4 1+(3N/4) 3N/4 
6 1+(5N/8) 5N/8 
8 1+ (35.N/64) 35N/64 
APPENDIX II 


Formal Solution of f(0) = /o"g(v)K(Ov)dv 


Here we shall find the formal solution g(v) of the integral 
equation 


f(0) = f “g(v)K (@v)dv, (A) 


when f(@) and K(6v) are assumed to be known functions. 
The Fourier transform method to be used is modeled on 
that of Payley and Wiener’® (the formal solution is also a 
simple consequence of some operations with Mellin 
transforms of f, g, and K." Let us make the substitutions 


6=e", yv=e (B) 
in (A) yielding 
ele) = f° glerer"K (e#")da. (C) 
Multiplying by 
ee "dn// 24 


and integrating: 


=f" e-"f (e-")e-iundn 


J/ 24J — «© 


=~}, ” g(e*)e~‘a"dax f * e—Ne-in(-@) K (¢—") dn, 
V £T J — oo —— 


Now let n—a=8 
1 © 
—n —1)\ p—tu 
Trad BR fle de ™ 
r Tyg fale tedaa fe Beh K (eA). 


10R. Payley and N. Wiener, Fourier Transforms (New 
York, 1934), p. 38. 

"Titchmarsh, Theory of Fourier Transforms (Oxford, 
1937), p. 315. 
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Define 
I(u) = fe PnK (e*)dp 
= f “xK (x)dx. 
0 
Then 
ee I(u) po ; 
cnieaiauas ” N\ piu gee exten —iau 
tm x. (ee a= Ne __ bee da 
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and taking Fourier transforms 


1 poe™edu po 
a =— pas 7 —@ “Hilti 
g(e*) = Toy JS Je \dn. 





By applying (B) 


1 po du po 
g(v) =5- Saw , £(8)(v)d8. 
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The Viscosity of Certain ‘Ferric Oxide’? Hydrosols *} 


Leon M. GREENSTEIN AND ARTHUR W. THOMAS 
Columbia University, New York, New York 
(Received December 11, 1941) 


(1) The changes produced in the viscosities of ‘ferric 
oxide” hydrosols by aging, heating, and the addition of a 
number of electrolytes and a non-electrolyte have been 
studied, in some instances together with the attendant 
changes in the pH values and specific conductances of the 
systems. (2) Agreement with Poiseuille’s law has been 
demonstrated for typical cationic and anionic hydrosols. 
(3) The decrease in viscosity of hydrosols on aging and 
heating is accompanied by an increase in specific con- 
ductance. This conductance increase in the case of cationic 
sols is due almost entirely to the increase in hydrogen ion 
activity which occurs during the reaction. This is not so 
for anionic hydrosols, where the conductance increase is 
explained by the conversion of coordinated citrato groups 
to citrate ions. (4) The addition of salts and acids to these 
hydrosols first results in a viscosity decrease, followed by 
an increase in viscosity at higher electrolyte concentra- 
tions. (5) The viscosity increase occurs at lower salt con- 
centrations for those salts which are potent precipitants 
than for those which are weak. For cationic hydrosols, salts 
of both weakly and strongly coordinating anions have the 
same effect in depressing the viscosity; in the case of 
anionic hydrosols, potassium citrate is less effective in 


INTRODUCTION 


HE addition of electrolytes to hydrosols 
has long been known to cause changes in 
the viscosities of the colloidal system. Farrow! 
showed that the viscosity of hydrosols of sodium 
palmitate is decreased by the addition of sodium 


* Presented at the Atlantic City meeting of the American 
Chemical Society, September, 1941, and published with 
the permission of the Society. 

} Dissertation submitted by Leon M. Greenstein in 
partial fulfillment of the requirements for the degree of 
Doctor of Philosophy in the Faculty of Pure Science, 
Columbia University. Publication assisted by the Ernest 
Kemptom Adams Fund for Physical Research of Columbia 

niversity. 


1F. D. Farrow, J. Chem. Soc. 101, 347 (1912). 


this respect than potassium sulfate and nitrate. (6) The 
effects of acids in decreasing the viscosity of a sol, either 
cationic or anionic, was found to be the same for all acids 
used, when compared at the same pH values. The more 
powerful precipitants caused the subsequent viscosity in- 
crease to occur at lower electrolyte concentrations. (7) 
Acids are less effective than salts in decreasing the vis- 
cosities of cationic hydrosols. (8) Large concentrations of 
acetone produced a decrease in the relative viscosity of 
cationic hydrosols; this effect is, however, negligible in 
comparison with the magnitude of the changes caused by 
electrolytes. (9) The viscosity decrease produced by aging, 
heating, and the addition of electrolytes may be attributed 
in each case partly to the increase in electrolyte concen- 
tration of the system. (10) The viscosity resulting from 
these reactions is considered to be dependent also upon the 
micellar charge: the lower the charge, the lower the vis- 
cosity. (11) The viscosity increase at higher electrolyte 
concentrations may be ascribed to agglomeration of dis- 
persed particles. (12) A mechanism based upon a Donnan 
equilibrium between the dispersed particle and the dis- 
persion medium has been proposed to account for the 
effects observed. 


chloride or sodium hydroxide, but passes through 
a minimum, increasing with further increase in 
electrolyte concentration. Woudstra? demon- 
strated that curves of similar shape were ob- 
tained upon the addition of electrolytes to ferric 
oxide hydrosols, and Pauli’ found electrolyte 
addition to lower the viscosity of ceric oxide 
hydrosols. In the ferric oxide system, viscosity 
decreases were obtained also on heating, and in 
the ceric, on aging. The extensive investigations 
of Loeb‘ on gelatin indicated that the depression 

2H. W. Woudstra, Kolloid Zeits. 8, 73 (1911). 

3 W. Pauli and A. Fernau, Kolloid Zeits. 20, 20 (1917). 


‘J. Loeb, Proteins and the Theory of Colloidal Behavior 
(McGraw-Hill Book Co., 1924), second edition. 
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of the viscosity of gelatin chloride dispersions 
was dependent upon the valence of the anion of 
the added salt: the greater the valence, the 
greater the depression. Similarly, the viscosity 
decrease for dispersions of gelatin on the alkaline 
side of the isoelectric point was found to depend 
on the valence of the cation of the added 
electrolyte. 

More recently, Thévenet and Boutaric have 
demonstrated the influence of the valence of the 
cation of the added salt on the lowering of the 
viscosity of arsenious sulfide hydrosols (anionic 
micelles) ,° while showing the valence of the anion 
to be of importance in the case of cationic 
micellar ferric oxide sols.* In the latter case, as 
in that of gelatin chloride,‘ the greater depression 
resulted upon the addition of the anion of higher 
valence. 

The problem of the viscosity of dispersions 
was early considered by Einstein, who derived 


- the equation 


nr=n/no=1+2.5¢ (1) 


for dilute dispersions of non-charged, very small, 
rigid spheres, which are large in comparison with 
the molecules of solvent.” * Here 7, is the relative 
viscosity, 7 is the viscosity of the dispersion, 
no the viscosity of the medium, and ¢ the ratio 
of the volume of the dispersed phase to the 
volume of the dispersion. The constant 2.5 
assumes greater values for particles which have 
shapes other than the spherical. 

Smoluchowski® presented an equation for the 
case of charged spheres by adding a term to 
correct for the ‘‘electroviscous’”’ effect : 


3 Diy? 
nr=1425¢| I+ (—) | (2) 
2rnor?\ 2x 


In addition to the symbols used above, dX is the 
specific conductance of the dispersion, 7 the 
radius of the spherical particle, D the dielectric 
constant of the medium, and ¢ the “electro- 


kinetic potential of the double layer.” The 
3 





factor $ of the final term did not appear in 


5A. Boutaric and S. Thévenet, Comptes rendus 205, 
862 (1937). 

6S. Thévenet, Comptes rendus 207, 68 (1938). 

7 A. Einstein, Ann. d. Physik (4) 19, 289 (1906). 

8 E. Guth, Kolloid Zeits. 74, 147 (1936). 

*M. V. Smoluchowski, Kolloid Zeits. 18, 190 (1916). 


AND A. W. 





THOMAS 


Smoluchowski’s original equation, but was added 
in a re-derivation by Krasny-Ergen.'® The 
conditions necessary for the derivation, as given 
by the latter, include: 

1. The tiny spheres must have uniformly 
charged surfaces. 

2. They must be complete insulators. 

3. The diffuse layer, or ionic atmosphere, must 
be very thin in comparison with the radius of 
the particle. 

The equation predicts that, for constant 
values of y, the viscosity of the dispersion will 
be greater the greater the charge on the particle 
and the smaller the conductance of the dispersion. 
The equation has been used by Kruyt and his 
co-workers" ” to explain the decrease in viscosity 
produced by added electrolyte, since this addi- 
tion both decreases ¢ and increases \. Considera- 
tion of these two variables alone, however, 
represents an oversimplification of the systems 
which they investigated, because these were 
hydrophilic, rather than hydrophobic, disper- 
sions. The particles are expected to be highly 
hydrated, and no assurance can be obtained that 
the hydration, and therefore g, will not vary in 
the presence of electrolytes. Kruyt and Tendeloo” 
themselves found that the viscosity and ¢ for 
gum arabic dispersions varied in opposite direc- 
tions, and assumed therefore that g was not 
constant. They suggest that the dispersed 
particles are surrounded by a shell of water, 
whose thickness changes in the presence of 
electrolytes. 

The Smoluchowski equation has been applied 
by Briggs'® to measurements on sodium gum 
arabic and caseinate dispersions, which contained 
various concentrations of sodium chloride. Values 
of no/(n—no) were plotted against ¢/[A(n— 7) ], 
since this choice of coordinates should yield, if 
hydration changes do not occur, a straight line 
whose slope is dependent only upon 7*. Although 
straight lines were obtained, the calculated 
value of the radius of the particle was found to 
decrease as the concentration of the dispersed 
phase increased. This variation in the value of 
the radius is not in agreement with the equation. 
10 W. Krasny-Ergen, Kolloid Zeits. 74, 172 (1936). 
1H. G. B. de Jong, Rec. trav. chim. 43, 189 (1924). 

”® Kruyt and Tendeloo, Kolloidchem. Beihefte 29, 413 


(1929). 
3D. R. Briggs, J. Chem Phys.. 45, 866, 943 (1941). 








VISCOSITY 


OF 


HYDROSOLS 


TABLE I. Description of hydrosols. 











Method 
of prep. 


Micellar 


charge Appearance of sol 


Time of meq 
dialysis* 


Total 

conc. conc. 
Fe: Anion 

meq. 

X-/l. 


Fet+t*/1. viscos.T 





Pos. Dark red-brown, 
slightly turbid to re- 
flected light 

Clear, very dark red- 
brown 

Similar to sol 1 
Clear, very dark red- 
brown 

Clear dark red-brown 
Clear red-brown 
Clear red-brown 
Very turbid red- 
brown by reflected, 
clear deep red by 
transmitted light 
Clear dark red-brown 
Clear dark red-brown 


Pos. 


Pos. 
Pos. 


Pos. 
Pos. 
Pos. 
Pos. 


Neg. 
IV Neg. 


8 days 


14 days 


16 days 
3 days 


5 days 
31 days 
31 days 
17 hours 
at 70°C, 
5 days at 
room temp. 
2 days 
36 hours 


1,219 54.9 1.570 


935 1.37 
556 
1,898 


546 
726 
641 
744 








* Approximately four liters of distilled water were used per day. 
tFresh sol, diluted: 1 ml! water to 10 ml sol. 


Loeb has made use of change in hydration as 
the mechanism for viscosity variations in gelatin 
dispersions. He determined that, using pH 
values as abscissae, curves of the same shape 
were obtained whether the ordinate was swelling 
of gelatin granules, or the osmotic pressure or 
viscosity of gelatin dispersions. All these proper- 
ties are at a minimum at the isoelectric point. 
The osmotic pressure curve was explained on the 
basis of a Donnan equilibrium, and the swelling 
and viscosity curves by considering the distribu- 
tion of ions which result from such an equilibrium 
when the charged gelatin particle is immersed 
in an electrolyte solution. The Donnan equi- 
librium predicts the hydration of the particle to 
be greatest when the charge is greatest and the 
electrolyte concentration of the solution smallest. 
Since the viscosity of the dispersion will vary in 
the same direction as the hydration of the 
dispersed phase, these are the same conditions 
as those required for maximum viscosity by the 
Smoluchowski equation. The mechanisms, how- 
ever, are quite different. 

The increase in viscosity which occurs when 
greater quantities of electrolyte are added is 
considered to be due to an agglomeration of 
individually dispersed particles. Smoluchowski 
points out that, for spherical particles, even the 
closest possible packing must produce an increase 


in the relative volume of the dispersed phase.° 
Should a looser packing occur, with more 
occluded water, the increase in relative volumes, 
and consequently in viscosity, must be still 
greater. Freundlich" has used this increase in 
viscosity as a means of measuring the rate of 
coagulation of aluminum hydroxide sols on the 
addition of various salts. 

Reports on the viscosity of “iron oxide” 
hydrosols*? *15!6 have been brief, and do not 
afford a knowledge of the factors upon which 
the viscosities of these systems depend. It is the 
purpose of this investigation to determine these 
factors by the study of a number of variables 
which influence the viscosities of these sols, and 
to propose a mechanism through which the 
changes in viscosity occur. 


EXPERIMENTAL 


Methods of Preparation of Hydrosols 


I. (A) One mole of ferric chloride was treated 
with excess ammonia. Each portion of the 
precipitate, after dividing among four half-liter 


14H. Freundlich and C. Ishizaki, Trans. Faraday Soc. 9, 
66 (1913). 

16 0. H. Ayres and C. H. Sorum, J. Phys. Chem. 34, 2826 
(1930). 

16 J. H. Yoe and E. Freyer, J. Phys. Chem. 30, 1389 
(1926). 
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Pressure, cm. water at 25°C 


Fic. 1. Variation of time of flow of sols 2 and 9 
with pressure. 


centrifuge cups, was washed centrifugally with 
eight portions of distilled water of 250 ml each. 
The precipitate was then suspended in distilled 
water, the volume of the suspension being one 
liter. Peptization was effected by adding either 
(a) 1M ferric chloride or (b) concentrated 
hydrochloric acid, to the extent of approximately 
0.1 equivalent of added electrolyte per equivalent 
of precipitated iron. (B) The same procedure 
was followed using ferric nitrate, and peptizing 
with 1M ferric nitrate. About 0.25 equivalent 
of ferric nitrate was added for each equivalent 
of precipitated iron. 

II. To an aqueous 2M solution of ferric 
chloride was added slowly, with agitation, 
concentrated ammonia equivalent to two-thirds 
of the ferric ion. 

III. 0.5 liter of 0.5M solution of ferric nitrate 
was added dropwise with agitation to one liter 
of distilled water at 70 to 80°C. Dialysis was 
carried on simultaneously to remove the products 
of hydrolysis. 

IV. A dialyzed sol (to absence of thiocyanate 
test for ferric ion on the diffusate) prepared by 
method II was precipitated and redispersed by 
adding, with agitation, potassium citrate solution. 

Dialysis in each case was accomplished by 
placing the sol in a nitrocellulose bag which was 
immersed in a beaker through which distilled 
water was allowed to flow continuously. 

The hydrosols are described in Table I. 
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Apparatus 
Viscosity 


A Bingham viscometer,’ having a working 
volume of five milliliters, was used for the 
determination of outflow times of hydrosols 
under various pressures. The efflux time for 
water at 25.000+0.005°C was 215.0 seconds 
under a pressure of 164.0 cm of water. For all 
other viscometric measurements, Ostwald vis- 
cometers were used. At 25°C, these had outflow 
times of about ninety seconds when containing 
five milliliters of water. All measurements of 
viscosity were carried out at 25°+0.005°C. 


pH Values 


These were measured at 25.0+0.1°C by the 
glass electrode. 


Conductance 


Twenty-five ml of the sol were placed in an 
unplatinized platinum electrode conductivity 











Spec. Conductance 
Mhos -/0# 
\% N 
& 3 
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N) 
8 
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felative Viscosity 
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4 ,-] 42 46 20 24 
Hours Heated at 80°C 


Fic. 2. Effect of heating on sol 6 (cationic). 


17E. C. Bingham, Fluidity and Plasticity (McGraw-Hill 
Book Co., 1922), first edition, p. 76. 
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cell. The cell constant was determined at various 
resistances by the use of several concentrations 
of potassium chloride, the equivalent conduct- 
ances of which are given by Shedlovsky."® 
Measurements were made at 25.000+0.005°C, 
using 1000-cycle alternating current and the 
modified Wheatstone bridge arrangement de- 
scribed by Dike.'® 


Determination of Nature of Flow 


According to Poiseuille, 
° n= Cpt—C'p/t, (3) 


where 7 is the viscosity of the liquid, p is the 
pressure causing flow, ¢ is the time of outflow of 
a given volume through the capillary, C and C’ 
are constants for a given viscometer, and p is 
the density of the liquid.*° In order to show 
whether the sols investigated obeyed Poiseuille’s 
law, their times of outflow through a Bingham 
viscometer were determined under various 
driving pressures. For large values of ¢, the 
second term becomes negligible, and 


p=k(1/t), (4) 


where k=7/C. The plot of 1/t¢ against p in 
Fig. 1 demonstrates that typical cationic and 
anionic sols* exhibit viscous flow. 



















4310 
4290}- 
2 
> 
Q 1.270h 
8 
g 
Lt 4250 F- 
; 
4230 |- 
OHected So/ 6. 
@Unhected So/ 6 plus 
hydrochloric acid. 
42/10 l l | i 
$0 4.6 4.2 J.8 34 FO 
pH 
Fic. 3. 


‘ST. Shedlovsky, J. Am. Chem. Soc. 54, 1411 (1932). 

' P. H. Dike, Rev. Sci. Inst. 2, 379 (1931). 

*° Reference 17, p. 74. 

*For brevity, hydrosols with positive micelles will be 
described simply as ‘‘cationic sols,”’ and those with negative 
micelles as ‘‘anionic sols.” 
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Fic. 4. Change in viscosity of sol 4 on aging at 25°C. 


TABLE II. pH value and specific conductance 
of heated sol 6.* 











1 2 3 4 5 6 7 
Hr. 
heated 
at 
80°C pH Qyz+°10® Aggy + 10° Akoaie 10° Kegic 10° Kobs *10° 
0.0 4.63 2.34 22.7 
0.5 4.42 3.80 1.46 0.51 23.2 23.0 
1.5 4.33 4.68 2.34 0.82 23.5 23.1 
3.5 4.19 6.46 4.12 1.44 24.1 23.6 
6.5 4.11 7.76 5.42 1.90 24.6 24.2 
11.5 4.00 10.00 7.66 2.68 25.4 25.5 
24.0 3.92 12.92 9.68 3.39 26.1 26.2 
72.0 3.70 20.0 17.6 6.17 28.9 29.0 





| 





« =specific conductance 
Ady + =4y + —2xy+,; ayz+, =2.34 10% 
Axcale = (0.350) Aayy+ 
Keale = +K0; ko =22.7 -1075 
* Sol diluted prior to heating: 10 ml water to 100 ml sol. 


Heating and Aging 


Portions of sol 6 (cationic), sealed in steamed 
seventy-five-ml Pyrex test tubes,t were placed 
in an 80+0.3°C-water thermostat. Test tubes 
were withdrawn from the bath at specific time 
intervals, and rapidly cooled in an ice-bath. 
The viscosity and pH value fall with time, 
while the conductance increases, all measure- 
ments being made at 25°C (Fig. 2). 

The lower curve of Fig. 3 shows this viscosity 
and pH value data plotted against one another. 
The addition of hydrochloric acid to the un- 
heated sol also lowers the viscosity and the pH 
value of the sol, and the effects of increasing 
concentrations of hydrochloric acid are shown by 
the upper curve of Fig. 3. The viscosity of the 
heated sol, at a given pH value, is seen to be 
distinctly lower than that of the unheated sol. 

It is interesting to note that the change in 
specific conductance on heating may be ac- 
counted for largely on the basis of the change in 


+ Samples heated in tubes which had been treated with 
sulfuric acid-dichromate cleaning mixture, but had not 
been subsequently steamed, gave erratic results. 
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Fic. 5. Aging of sol 9 (anionic) at 25°C. Lower pair 
of curves: the aging at 25°C of a portion of the sol which 
had been heated at 90°C for 13.5 hours. Solid lines: 
viscosity. Broken lines: resistance. 


TABLE III. Sol 4 (undiluted). 


Age in days* 0 315 
Rel. viscosity 1.843 1.587 
pH 3.46 3.17 


8.92 x 10-4 10.35 X 1074 


* Refers to time elapsed after removal from dialyzing bag. 


Spec. conduct. (Mho) 


hydrogen ion activity. Column 4 of Table II 
shows the increase of hydrogen ion activity of 
the heated portions of sol over the hydrogen ion 
activity of the original unheated sol, as calculated 
from the pH values. This increase in hydrogen 
ion activity should lead to the increase in 
specific conductance shown in column 5, as- 
suming 350, the equivalent conductance of the 
hydrogen ion at infinite dilution (25°C),”! to be 
the equivalent conductance of the hydrogen ion 
in this system. In column 6, these conductance 
increases are added to the specific conductance 
of the original sol, and in the last column are 
given the observed conductances for the heated 
portions of the sol. 

The viscosity of cationic sols aging at room 
temperature follows a path similar to the heating 
curve discussed above. Figure 4 shows the 
change of viscosity of sol 4 with time.t That 
the pH value and specific conductance of the sol 
also change is shown by Table III. 

21 Landolt-Bornstein Tabellen (Julius Springer, 1936), 
Vol. 3, p. 2059. 


t Sol diluted just before making measurements: 10 ml 
water to 100 ml sol. 
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The decrease, with time, of the viscosity of an 
anionic sol is shown in Fig. 5.{{ Heating a 
portion of this sol causes a much more rapid 
decrease in viscosity: 13.5 hours at 90°C causes 
the relative viscosity of the sol to fall from 1.273 
to 1.212, and the pH value from 6.67 to 5.97, 
while the specific conductance increased from 
5.25-10- to 6.90- 10-4 mho, a change of 1.65-10-* 
mho. The change in hydrogen ion activity 
(8.6-10-7 mole/liter) will account for a change 
in conductance of only 0.0030-10-* mho. Thus 
it is necessary to conclude that ions other than 
hydrogen are mainly produced. These must be 
citrate ions. 

It is to be noted that the aging of the heated 
portion of sol 9 produces some reversal of the 
effect of the heating, Fig. 5 indicating a small 
increase in viscosity and in electrical resistance. 


Effect of Added Salts 


To ten ml of the sol in a small glass-stoppered 
bottle was added one ml of salt solution. In 
Table IV are given the viscosities resulting from 
the addition of several concentrations of potas- 
sium sulfate solution to sol 4. Column 1 lists 
the viscosity immediately after the addition of 
the salt solution, and columns 2 and 3 the 
viscosities at subsequent intervals. 

The aging of the sol to which potassium 
sulfate has been added is accompanied by a 
viscosity decrease. However, the decrease in 
viscosity with time is smaller for the sol—salt 
mixture than for the sol diluted only with 
water. 

If larger concentrations of salt are present, or 
if the sol is one whose viscosity does not diminish 
rapidly with time, another effect becomes 
apparent. Figure 6 shows the aging of sol 3, in 


TABLE IV. Addition of potassium sulfate to sol 4. 








Relative viscosity 





Immediately 72 hours 216 hours 
Final conc. of after after after 
pot. sulfate addition addition addition od 
0.00 meq./1 1.731 1.703 1.663 
0.91 1.600 1.575 1.549 
1.82 1.535 1.518 1.495 
3.64 1.442 1.438 1.432 











tt Aging experiment conducted on diluted sol: 10 ml 
water to 100 ml sol. 
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the presence of added potassium sulfate and 
potassium chloride. The viscosity increases with 
time, and the rate of increase is greater where 
the salt concentration is greater (curves B and 
C). It is evident also that for equal salt concen- 
trations (curves C and D), the rate of viscosity 
increase is greater for potassium sulfate than for 
potassium chloride. 

In order to obtain comparable viscosity values 
for viscosity-salt concentration curves, all values 
for a given experiment were determined either 
immediately after salt addition, or twenty-four 
hours after salt addition.* The viscosities in the 
presence of a given concentration of all salts 
were determined within a time sufficiently short 
to render insignificant differences due to the 
aging of the sol. 

The effect of salts on the viscosity of cationic 
sols is shown in Figs. 7 and 8, and on anionic 
sols in Fig. 9. 

In these figures, the upward trend in the 
curves is attained at lower salt concentrations 
for salts which are potent precipitants than for 
those which are weaker in this respect. (For 
sols in which the dispersed phase is positively 
charged, the potent precipitants are the salts of 
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Fic. 6. Change of viscosity of sol 3 (cationic) with time 
at 25°C. A. No salt added. B. Added potassium sulfate, 
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. Data for sols 8 and 9 were obtained immediately after 
salt addition; data for sols 1, 5, and 10 were obtained 
twenty-four hours after salt addition. 
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Fic. 7. Effect of potassium salts on viscosity of sol 1 
(cationic). 


anions which are good coordinators with iron. 
For anionic sols, the order of precipitating 
ability is reversed, since strongly coordinating 
anions may augment the negative charge of the 
micelle, thus conferring greater stability upon 
the sol.) The decrease in viscosity, however, 
appears to be independent of the anion used in 
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Fic. 8. Effect of potassium salts on viscosity of cationic 
sols 8 (upper) and 5 (lower). Points corresponding to 
viscosities outside the range of the plot are shown by 
arrows. The abscissa is indicated by the position on the 
graph, and the ordinate by the number beneath the arrow. 
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the case of positive sols. The curves shown in 
Figs. 7 and 8 are representative of those obtained 
for all positive sols in this research, whether 
preparation was by method I, II, or III. These 
disagree with the data of Thévenet,® who 
reported the decrease in viscosity to be in the 
order K,Fe(CN).> K3PO, > K.SO,; > KCl. There 
is evident, nevertheless, a difference in the effect 
of various anions on the decrease of viscosity of 
anionic sols. Potassium citrate produces a smaller 
viscosity decrease than salts of the less strongly 
coordinating sulfate and nitrate ions. 


. Effect of Added Acids 


The above procedure was followed, except 
that acid additions replaced the addition of salts. 
Figure 10A shows the results on plotting, as 
above, relative viscosity against final electrolyte 
concentration. (The acid curve which is included 
in Fig. 7 shows the relative efficiency of acids 
and salts in lowering the viscosity.) The order of 
the curves of Fig. 10A indicates that the stronger 
the acid, the greater the viscosity depression, 
which suggests that the smaller depressions are 
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Fic. 9. Effect of potassium salts on viscosity of anionic 
sols 9 (lower) and 10 (upper). 
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due to smaller actual ion concentrations. This 
was tested by plotting the viscosity of each sol- 
acid mixture against its pH value (Fig. 10B). 
At constant pH value, the viscosity decrease for 
this positive sol is apparently independent of the 
nature of the anion. Again, the strongly co- 
ordinating anions cause the viscosity increase to 
appear at lower electrolyte concentration. Similar 
data were obtained with an anionic sol (Fig. 11), 
there being a difference owing to various anions 
only for the ascending portions of the curves. 
The viscosity minimum is first produced by the 
acid with the most weakly coordinating anion, 
as in the case of salts on this sol (Fig. 9). A 
comparison of the effects of an acid and a 
potassium salt on sol 10 may be seen in the 
nitrate curves of Fig. 9. 


Effect of Neutral Molecules 


To forty ml of sol were added ten to sixty ml 
of acetone, and the volume made up to 100 ml. 
The relative viscosity of each solution was 
obtained by dividing the viscosity of the sol- 
acetone-water mixture by the viscosity of the 
acetone-water solution of the same acetone 
concentration. Although the viscosity is lowered 
by the addition of acetone (Table V), the effect 
is negligible when compared to the viscosity 
decrease produced by an electrolyte. The 
relative viscosity of this sol is reduced, for 


TABLE V. The effect of acetone on sol 7. 











Conc. acetone Rel. 
(moles/1) viscosity 
0.00 1.144 
1.35 1.141 
2.70 1.137 
4.05 1.134 
5.40 1.130 
6.75 1.127 
8.10 1.123 








example, from 1.144 to 1.031 in the presence of 
0.004M potassium nitrate. 


DISCUSSION 


According to the mechanism proposed by the 
workers of this laboratory,” one of the types of 


22 A. W. Thomas and co-workers, J. Am. Chem. Soc. 60, 
2384 (1938) ; 58, 2526 (1936) ; 57, 1825, 2131 (1935) ; 57, 1821, 
2538 (1935); 57, 4 (1935); J. Phys. Chem. 35, 27 (1931); 
Kolloid Zeits. 86, 279 (1939). 
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chemical changes occurring during the aging or heating of metal oxychloride dispersions may be 


formulated :t 
H.0 Cl OH HO H,0]** 
we Dae 
Fe Fe 
a Py 


H.O H.0 OH’ HO “H.O 


This reaction explains the increase in the hy- 
drogen ion activity and in the conductance of the 
sol. The experiments discussed above indicate 
that the addition of any electrolyte to an iron 
oxychloride, -nitrate, or -citrate sol diminishes 
the viscosity of the dispersion. The viscosity de- 
crease which accompanies the heating and aging 
reaction may thus be accounted for partly by the 
increased hydrochloric acid concentration of the 
sol. A simultaneous result of the above reaction is 
a decrease in particle charge. Increased electro- 
lyte concentration may be realized without the 
decrease in charge by the addition of hydrochloric 
acid to the unheated sol. The micellar charge may 
in fact be increased by reversing the above reac- 
tion. (Hazel and Ayres” have shown that the 
electrophoretic velocity of the particles in posi- 
tive ferricoxychloride sols increases with the de- 
crease in pH value produced by the addition of 


* The micelles of metal oxide dispersions are considered 
to be complex coordination structures containing a large 
number of metal atoms, linked to one another by means 
of ol or oxo bridges. Other coordination positions are 
occupied by anions, and by aquo and hydroxo groups. 

alt Hazel and G. H. Ayres, J. Phys. Chem. 35, 2930 


H:O Cl OH HO OH ]** 


+ H.O—H;0 i Fe Fe 


/ | J | 
H,O H»O OH H:2O H.O]. 


acid, suggesting an augmented positive charge.) 
At equal pH values, the heated and unheated sols 
have the same hydrochloric acid activities, and 
differences in the viscosity should be attributable 
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to differences in micellar charge. Figure 3 
demonstrates that the heated sol has, at a 
given pH value, a lower viscosity than the 
unheated, indicating that the viscosity is 
decreased when the micellar charge is lowered. 

Another illustration of this point may be seen 
in Fig. 7, where the viscosity diminution pro- 
duced by nitric acid is less than that caused by 
potassium nitrate, because of the increase in 
particle charge by reaction of the hydrogen ion 
with the micelle. A contributing factor in this 
instance, however, is the decrease in electrolyte 
concentration of the dispersion which results 
from the combination of some of the added 
hydrogen ion with the micelle. 

The curves for the aging of the negative 
dispersion (Fig. 5) may be considered to be 
analogous to those for the positive, the main 
reaction being the hydrolysis of the micelle in 
such a manner as to replace coordinated citrato 
groups with aquo groups. That this occurs was 
indicated by the change in conductance which 
occurs on aging. This reaction decreases the 
negative charge of the particle and increases 
electrolyte concentration, producing the decrease 
in viscosity. The increase in viscosity and 
electrical resistance on aging after heating may 
be the result of recoordination of some of the 
citrate ions at the lower temperature. 

Since the dependence of the viscosity of 
“ferric oxide” hydrosols on charge and electrolyte 
concentration has now been demonstrated, 
certain effects must be expected on the addition 
of salts to the dispersion. All salts will increase 
the electrolyte concentration; thus all salts may 
be expected to decrease the viscosity. Those 
salts which have strongly coordinating anions 
must decrease the charge of positive micelles to 
a greater extent than those with poorly coordi- 
nating anions, while those same salts will tend 
to augment the charge of negative micelles more 
than will the salts of the poorly coordinating 
anions. In the case of positive micelles, therefore, 
the better coordinators should produce a greater, 
and in the case of negative particles, a smaller, 
decrease in viscosity. The latter expectation is 
fulfilled (Fig. 9), while the former, in the present 
research, is not (Figs. 7 and 8). It will be seen 
in all these curves that the minimum in the 
viscosity-electrolyte concentration curve is at- 


AND A. W. 





THOMAS 








tained at lower electrolyte concentrations for 
potent precipitants than for poor. Herein lies 
an explanation for the apparent similarity of 
effect of all anions on the viscosity decrease of 
cationic sols. The strongly coordinating anion 
reduces the positive charge of the micelle. This 
should produce a decrease in viscosity. The 
stability of the sol has at the same time been 
decreased by the decrease in charge, with the 
result that agglomeration of particles may occur, 
thus increasing the viscosity. (It was, in fact, 
possible in some instances to detect by eye very 
minute suspended particles in the citrate series, 
as the minimum in the curve was approached.) 
An indication of the effect of agglomeration may 
be seen in the slow increase of viscosity with 
time in Fig. 6. Where the salt concentrations 
are equal, as in curves C and D, the more 
strongly coordinating anion causes the more 
rapid viscosity increase. This balance of opposite 
forces produces a smaller decrease in viscosity 
upon the addition of good coordinators than 
would be the case were no agglomerating action 
present. Since both the viscosity-decreasing and 
the viscosity-increasing reactions have the same 
source (the lowering of the micellar charge by 
the coordination of anions), the variation with 
extent of coordination does not appear, and all 
the salt curves tend to coincide. 

The same argument may be applied to the 
acid curves. Differences among the effects of 
various acids are obvious only at the ascending 
portions, and it is the tendency toward agglom- 
eration in the presence of potent precipitants 
that masks the viscosity-decreasing effect at the 
lower acid concentrations. 

In the one instance where differences in the 
descending parts of the curves are produced by 
different anions (the addition of potassium salts 
to anionic sols, Fig. 9), an electrolyte concen- 
tration factor, as well as a particle charge factor, 
may contribute to the result. The coordination 
of citrate ions will lower the actual electrolyte 
concentration of the dispersion, and thus the 
final equilibrium potassium citrate concentration 
will be smaller than the final salt concentration 
in the case of an equal quantity of added nitrate. 

That the viscosity changes of the magnitude 
observed are due only to electrolyte influence is 
shown by Table V above, where acetone up to 
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sixty percent by volume caused a relatively 
small change in the viscosity. This viscosity 
decrease may possibly result from a lowering of 
the particle charge by increased association (in 
the Bjerrum sense*) of positive micelle and free 
chloride ions, through reduction of the dielectric 
constant of the dispersion medium by the 
addition of acetone. 


Systems described by the Smoluchowski Eq. 
(2) would show a dependence of viscosity on 
electrolyte concentration and micellar charge, 
through the presence of the variables \ and ¢. 
Freundlich®> points out, however, that iron oxide 
sols, in regard to their viscosity, are more closely 
related to the hydrophilic sols than to the 
hydrophobic. Since the hydration of the particles 
is great and unknown, no estimation of the value 
of g¢ of the Einstein (1) and Smoluchowski 
equations may be obtained, and the applicability 
of the equations to these systems cannot be 
ascertained. It is questionable, indeed, whether 
a ferric oxide dispersion is one composed of the 
rigid, non-conducting spheres, surrounded by 
relatively thin ionic atmospheres, for which the 
Smoluchowski equation was derived. According 
to the model to be presented below, which is 
based upon the properties of “ferric oxide” 
hydrosols, the particles would be neither rigid 
nor non-conducting. There has been evidence 
that the dispersed particles in these systems are 
non-spherical,?® and may be disks or platelets in 
aged iron oxide sols.2?7 Furthermore, the radius 
of the particle is believed not to be great in 
comparison with the thickness of its ionic 
atmosphere, in colloidal dispersions of low ionic 
concentration.?* Freundlich states that the water 
content of the particles and the resultant ‘“‘active 
volume”’ are nevertheless of great influence on 
the viscosity of sols, generally in the sense that 
with increasing water content the viscosity 
increases. It is therefore of interest to consider 
whether the extent of hydration may be related 


**R. M. Fuoss, Chem. Rev. 17, 27 (1935). 

*°H. Freundlich, Colloid and Capillary Chemistry 
(Methuen and Co., Ltd., 1926), trans. from the third 
German edition, p. 368. 

*° B. Lange, Zeits. f. physik. Chemie 132, 1 (1928). 

“7H. Diesselhorst and H. Freundlich, Physik. Zeits. 17, 
117 (1916). 

**D. C. Henry, Proc. Roy. Soc. A133, 106 (1931). 


to the electrolyte concentration of the dispersion 
and the particle charge. 

Let us assume a model for the positively 
charged micelle in a ferricoxychloride hydrosol. 
Asa result of the positive charge, the distribution 
of ions in the solution will be such that an 
atmosphere of predominantly negative ions will 
surround the positive micelle. The micelle con- 
tains water, and water is free to diffuse into and 
out of it. Browne”® has shown that all the water 
in a ferricoxychloride system acts as solvent 
water for added dextrose, as determined from 
the freezing point depression. Thus the ions of 
the intermicellar liquid must be considered free 
to diffuse into the water within the particle, and 
the ionic atmosphere, instead of originating at 
the particle boundary, is continuous from the 
intermicellar liquid through the particle. The 
larger the positive charge for a particle of given 
size, the more dense will this negative ionic 
atmosphere be, and the greater will be the 
concentration of negative ions in the liquid 
within and around the particle. 

In the system under discussion, the negative 
ions are chloride ions, and the concentration of 
these ions within, and in the immediate vicinity 
of, the micelle is greater than in the bulk of the 
dispersion medium. In addition to chloride ions, 
the dispersion medium contains hydrogen ions. 
At equilibrium, as has been shown by Donnan,*° 
the energy required to transport reversibly and 
isothermally 6n moles of H+ from a part of the 
intermicellar liquid (region 1) to a part of the 
dispersion medium which is within the particle 
(region 2) is equal to the energy which can be 
gained from the corresponding reversible and 
isothermal transport of 6” moles of Cl-. The 
diminution of free energy for this process is 
zero, hence 


ay*, 


6nRT In 


acl» 
4+énRT In——=0, (5) 


an*; acr; 


Qn*, Aci"; = an*, Acti“. 


where the subscripts denote the aforementioned 
regions. For simplicity, the equation of products 
may be approximated using concentrations in 


*F. L. Browne, J. Am. Chem. Soc. 45, 297 (1923). 
30 F, G. Donnan, Zeits. f. Elektrochem. 17, 572 (1911). 
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place of activities: whereas the decrease of the charge to zero gives 
(H+}(Cr-]=[H+}[Cr}. (7) P =[(4x?)!—2eJRT=0. (12) 






Now, if x=[H+]i=(Cl-]:, y=[H*].=the con- ee 
centration of chloride ions which balances these 
hydrogen ions, and z=the concentration of 
excess chloride ions in region 2 which are held 





Adding hydrochloric acid to a ferricoxychloride 
sol may produce some increase in 2, but x will 
increase to a greater extent, and once more 









by the positive charge of the micelle, [CI~]s iP. =0. 
=y+z, and ion 
x?=y(y+sz), (8) Achange in osmotic pressure will be followed by 






~2+(22-++4x2)3 the diffusion of water into or out of the particle, 
' (9) causing a change in particle volume and in the 
2 fraction of the total volume occupied by the dis- 
persed phase. Hence, a viscosity change occurs. 
A similar treatment for the addition of an 
electrolyte without common ions, such as 
potassium nitrate, also predicts a viscosity 
P= (2y+2—2x)RT. (10) decrease. This mechanism is consistent with the 
influence of micellar charge, since (11) indicates 
that the greater degree of swelling occurs when 
, z has the greater value. Also explained is the 
P=[(e+4x!)'— 2x RT. (11) lack of a primary effect on the addition of a 
From (11) it may be seen that when zis not zero non-electrolyte, acetone, which does not directly 
(this will be the case when the micelle is charged), vary the ion concentrations, but merely substi- 
the osmotic pressure will have positive values, tutes throughout the system for water. 








7 





The difference in osmotic pressure between the 
two regions under discussion depends on the 
difference in ionic concentrations, or 









Substituting in this expression the value of y 
from (9) gives 
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The Approximate Solution of Schrédinger Equations by a Least Squares Method 
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A method of approximation to the solution of a Schrédinger equation has been developed in 
which variation functions are used but no integrations are involved. The procedure involves the 
evaluation of the energy for a set of representative points in configuration space. The parameters 
in the variation function are then chosen by applying the condition that the mean square devi- 
ation of the energy from the average should be a minimum. — 







INTRODUCTION tions have been applied to more complicated mole- 
HE calculation of the electronic energy cules. It would be very desirable to have methods 
levels of molecules requires the use of that could be used successfully in such cases. 
approximation methods. The hydrogen molecule Of the many approximation methods used in 
ion, H.+, and the hydrogen molecule, Hz, are quantum mechanics? the best known are the 
examples of molecules that have been solved with =————— 
2 See, for example, L. Pauling and E. B. Wilson, Jntro- 


high accuracy.’ Only rather crude approxima- duction to Quantum Mechanics (McGraw-Hill Book Com- 
ae eaa aad pany, New York, 1935), Chaps. 6 and 7. An example of a 

1 For Het, e.g. the work of G. Jaffe, Zeits. f. Physik 87, recently developed method is the numerical integration 
535 (1934); for Hz, H. M. James and A. S. Coolidge, J. | method of G. E. Kimball and G. H. Shortley, Phys. Rev. 
Chem. Phys. 3, 129 (1935). 45, 815 (1934). 
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perturbation methods and the variation method. 
Perturbation methods are unsatisfactory for the 
calculation of electronic energy levels of mole- 
cules because the perturbations involved are 
much too large. On the other hand the variation 
method if pushed far enough will give results of 
sufficient accuracy. It is this method that was 
used by James and Coolidge! in their calculation 
of the hydrogen molecule. 

The application of the variation method to 
more complicated systems is both difficult and 
tedious. Good variation functions may easily be 
chosen, but there is no straightforward method of 
evaluating the required integrals. If simpler 
functions are used so that the integrals become 
reasonable, then the determination of the best 
values of the parameters becomes a lengthy 
process. In this paper is described an attempt to 
develop a method which makes use of a variation 
function but in which no integrations are re- 
quired. Numerical calculation here is more 
extensive than in the variation method but is 
put in such a form that machine methods may be 
used.* 


GENERAL FORMULATION OF THE METHOD 


Let the Schrédinger equation for the stationary 
states of the system be Hy= Ey where H is the 
Hamiltonian operator and y is a function of the 
coordinates. E is the energy parameter. Dividing 
both sides of the equation by y, 


Hy /p=E. (1) 


Since E is a constant, independent of coor- 
dinates, then if y is a solution, Hy divided by y 
must be a constant independent of coordinates 
even though Hy and y are both functions of the 
coordinates. 

Suppose that ¢ is an approximate solution. 
Hy/y will in general be not a constant but a 
function of coordinates. As ¢ is varied to become 
more nearly a solution Hg/g becomes more 
nearly a constant which is equal to the energy 
of the particular state represented by ¢. 

Define 

e=Ho/¢. (2) 

* Punched-card machine methods are labor-saving where 
a long calculation involves many similar steps. For 
applications of punched-card machines see W. J. Eckert, 


Punched Card Methods in Scientific Computation (Columbia 
University, 1940). 


e may be evaluated at any point without dif- 
ficulty because H is a differential operator. Now 
take as a criterion of the excellence of the ap- 
proximate function g the mean square deviation 
of « from its mean value. To find the best solu- 
tion, vary the parameters in ¢ until the mean 
square deviation of « from its mean value is a 
minimum. Using the subscript «4 to indicate the 
taking of mean values this condition may be 
written as 

Minimize ((€— €s)?) as (3) 


or after expansion 
Minimize ((€)m— €’). (4) 


The method of taking the mean value is not 
unique. For example ¢, may be taken over a 
finite set of points so that ew=)>>; giei/>: 2: 
where g; is a weighting factor, or « may be 
taken as an integral over all points in con- 
figuration space. A special case of €4 represented 
as an integral is where ¢*¢ is used as a weighting 
factor. Then 


a= f eteedr/ f oredr 


or introducing the definition of « 


ax f ottteds / f oredr. 


In this case €, is the usual quantum-mechanical 
expectation value or mean value of the energy 
for the system in the state yg. Where other 
weighting factors are used or where summation 
over a finite set of points is resorted to €y may 
still be interpreted as a mean value of the energy 
and « itself as the energy at a given point. As a 
result of this the condition (3) or (4) for the best 
function ¢ may be interpreted as requiring the 
energy to be as nearly as possible independent of 
the coordinates. 

In the limit, for a correct solution, ¢€ is a 
constant, and regardless of what weighting pro- 
cedure is used the mean square deviation will be 
zero. This fact leads to the conclusion that a 
satisfactory method for taking the mean would 
be to use the mean over a set of representative 
points rather than to integrate over all possible 
points. Naturally the number and distribution 
of points would determine the accuracy with 
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which the approximation is made. However, 
the use of a set of points replaces integration 
by summation and therefore may eliminate an ex- 
tremely difficult integration. It should be noted 
that this replacement of integration by sum- 
mation is not just a numerical integration. 
In choosing the points for the summation there 
is no need to spread them evenly or by suitable 
weighting to cause the summation to approach 
to an integration. 


THE METHOD APPLIED TO A LINEAR 
VARIATION FUNCTION 


g= Lm CrnPms 


where the C,,’s are coefficients in a linear com- 
bination of ¢,,’s. Since it is the energies of sta- 
tionary states that are to be calculated the C,,’s 
and g's can be restricted to real numbers and 
real functions respectively. Substitution in (2) 


yields 
ieee : CmlT om/ > m CnOm 
€(t) = Dim CulL gm(t)/ Lim Cuem(t), 


where the index 7 is used to refer to values of 
functions at a particular point and /7¢,,(7) means 
the value of Hg at point 7. 

In order to obtain the least square deviation 
it is convenient to calculate (e?),, and €,? as 
used in condition (4). The weighting factor will 
be taken to be p;¢(z) since this results in a great 
simplification of the formulas due to the possi- 
bility of changing the order of summation. Then 


_ Di pie? (ie 

— Dep) 
_ Xi bi(Lm Cnenli)) (Ln Col on(i)) 
Le Pi(Lm Cugm(6)) (En Cnen(i)) 


and interchanging the order of summation over 
4 and m and n, which is permissible since these 
are all finite summations, 


= y - CLs di Divm(t)H¢en(t) 
Xm Xn CnCn Xi Pigm(i) onli) 


Likewise, 
(e?)w= Doi Pig*(t)eP?/DX: pig?(i) 
P. Dem Lon CmCn Doi PilT em(i)H ¢n(t) 
Lm En CaCn Li Pigm(den(i) 


Let 


Env 








Eay 
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To simplify the expressions three definitions 
are introduced: 


Ls PiGm(t) ¢n(t) =Da, x; 
Di Pign(t)H¢, (i) = n,n. 
i pi on(t)He,(t) =Kn, ne 


Formulas (5) and (6) then become 
Cu = Dom p Chua, at hae as Culbaln 00 (8) 
(ey = om a Chilean al Ze Zs C0. Daw (9) 


Having chosen a set of functions ¢,, and a set 
of points 7 the quantities Dm,n, Hm,n, and Km,» 
can be calculated. There remains the deter- 
mination of the C,,’s such that Eq. (4) with the 
aid of Eqs. (8) and (9) is satisfied. In general 
several choices of the relative values of the C,,’s 
are possible corresponding to different energy 
states of the system. 


THE DETERMINATION OF THE LEAST 
SQUARE DEVIATION 


Minimizing the mean square deviation of «¢ in 
this method corresponds to minimizing the 
energy integral in the variation method. Whereas 
in the latter process the use of a linear variation 
function leads to a secular equation to be solved, 
the present method, if a similar procedure were 
followed, would lead to equations of great com- 
plexity. Instead, a method of successive approxi- 
mations may be used for the determination of 
the values of the C,,’s which minimize the mean 
square deviation of e. 

Given a set of approximate C,,’s, the corre- 
sponding éy and (e*),, and mean square deviation 
(€*)y — €y2 may be calculated from Eqs. (8) and 
(9). A single coefficient C, is then changed by 
slight positive and negative values and (€*) a — én” 
determined for each. A new value of C;, is then 
estimated which will roughly minimize (¢*)s— €x”. 
The procedure is repeated for each C,, in turn 
following through several cycles of all the C,,’s 
until (¢?)—€? shows no appreciable further 
decrease. The value of ey for this final set of 
C,’s is then the desired approximate energy. 
Computational details of this procedure are to 
be found in the appendix. 
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CONVERGENCE OF THE METHOD AND THE 
CHOICE OF POINTS 


That the convergence of the method is inti- 
mately bound up with the choice of points is 
shown by the following argument. For any 
arbitrary approximation function g, e will be a 
definite function of the coordinates. Now choose 
all the sets of coordinates or points for which e 
has some constant value E’. The values of ¢ 
associated with this set of points are all equal 
and the mean square deviation of ¢ is zero. 
Therefore, if the set of points had been chosen 
originally as a basis for the present method of 
approximation, the mean square deviation of « 
could be minimized to zero without the function 
y being equal to the true solution y and without 
the value E’ being equal to the correct value E, 
all of this despite the fact that the set of points 
may have been infinite in number. The method 
does not necessarily converge to the correct 
solution merely by taking a sufficiently large 
number of points. It is probable that for the 
process to converge to the correct solution the 
points must be so chosen that in the limit all 
possible points are included. It is obvious for 
this limiting case of all points included that the 
mean square deviation of ¢ can be zero only when 
yg is the correct solution and e=E. For suppose 
that the mean square deviation of ¢ is zero. Then 
e=Hgy/g must be constant, say E’. Therefore, 
Hy=E’g and ¢ must be a solution y and E’ a 
value of the parameter E. 

The best choice of points to be used could 
_ possibly be found by a method similar to the 
method of Gauss for the case of ordinary nu- 
merical integration.t However, physical intuition 
offers a simpler method. The form of the wave 
function, as long as it is constrained to satisfy 
the boundary conditions, is most important 
where its square is largest, i.e., where the 
probability amplitude is greatest. This, of course, 
will be the configurations where the electrons 
are not too far from the nuclei. However, points 
in configuration space corresponding to actual 
Superposition of electrons on nuclei are to be 
avoided because at such points for any ordinary 
approximation functions ¢ is likely to become 


tJ. B. Scarborough, Numerical Mathematical Analysis, 
p. ane wish to thank the reviewer for reference to Gauss’ 
method. 
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infinite or many-valued. In the practical appli- 
cation of the method probably the choice of 
points for a given problem will have to be deter- 
mined through experience with similar problems. 


APPLICATION OF THE METHOD TO THE 
HYDROGEN MOLECULE-ION 


The theory presented in the preceding sections 
gives no indication of the accuracy with which a 
problem may be solved after a set of points has 
been chosen. In this section is described a prac- 
tical application which was carried out for the 
purpose of determining the approximate number 
of points and functions necessary to arrive at a 
solution of sufficient accuracy for a problem of 
the type under consideration. The problem of 
the normal state of the hydrogen molecule-ion 
was chosen for this purpose because of its sim- 
plicity and because of the fact that an accurate 
solution is known. 

The Hamiltonian operator expressed in atomic 
units‘ is 


H= —(1/2)¥?—(1/ra) —(1/re), 


ra and rg being the distances of the electron 
from nucleus A and nucleus B, respectively. For 
the normal state y may be expressed as a function 
of r4 and rg as coordinates inasmuch as there is 
cylindrical symmetry about the internuclear 
axis. The operator V? may be expressed in terms 
of the coordinates 74 and rg together with a 
third, angular, coordinate determining position 
about the axis. All derivatives with respect to 
this angular coordinate vanish since y is a 
constant so far as this coordinate is concerned. 
In these coordinates the Schrédinger equation 
becomes 


1 ay 
2 Ora? 


1 d*y 
2 Orp? 
1 dy 1 od 1 
-——--——-— ¥-—v 


YA Ora ’B Orr Tr 


where cos AB is the cosine of the angle between 


4 Unit of length—radius of ist Bohr orbit of H; unit of 
mass—mass of electron; unit of time—(1/27) period of 
electron in first Bohr orbit of H. These result in the unit 
of charge being the electronic charge and the unit of 
energy being twice the ionization energy of an H atom 
with infinite nuclear mass. 
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¢ to 20 30 


Fic. 1. Ordinate: (€*)4,—e?, atomic units; abscissa: 
number of variations; dotted line: 6-point calculation; 
dashed line: 9-point calculation; full line: 15-point 
calculation. 


the radius vectors from the electron to the two 
nuclei. This comes in through the relation 


cos AB=(r4?+rp?—R?)/2rarz, 


where R is the internuclear distance. The value 
of R will be taken as 2, corresponding to the 
minimum in the potential energy curve. 

Although the equation is separable in elliptical 
coordinates and not separable in the coordinates 
ra and rg, the latter set will be used so as to have 
a system of calculation which can easily be 
generalized to more complicated molecules. 

y will be approximated to by the function ¢ 
of the form 


g=exp [—2(ratrs) ]F(ra, rp), 
F(ra, 1B) = Dee Doi Ceita* rp’. 
Substitution in Eq. (2) yields 


where 


—eF=>) & ratra!| 2 +cos AB) 


k,l 


1 1 z 
—(z— »(—+—) —k—(1+cos AB) 


Ya TB TA 


z 1 
—l—(1+cos AB)+k(k+1)—_ 


’R 2ra? 


1 cos AB 
+10+1)—+e———. (10) 


2rpz? TATB 


Since the lowest state of H:+ is symmetrical 
with respect to an interchange of r4 and rg, the 
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Fic. 2. Ordinate: — ey, atomic units; abscissa: number 
of variations; dotted line: 6-point calculation; dashed 
line: 9-point calculation; full line: 15-point calculation; 
— line at —es,=1.1026 indicates known correct 
value. 


set of approximation functions may be grouped 
as follows: 


gi=exp [—2(rat+rz) ], 

g2=exp [—2(rat+re) |(rat+rs), 

¢gs=exp [—2(ratre) ](ra2+rp?), at) 
gs=exp [—2(ratrs) Jrare, . 
gs=exp [ —2(ratra) |(ra?+rp°), 

gs=exp [—2(ratre) |\(ra*ret+rars’), 


etc. 


The set of six functions above are used in this 
calculation together with F? as the weighting 
factor. 

The parameter z could be varied as well as the 
coefficients Cm. However, in this calculation it 
was given a single constant value. A reasonable 
value of z was determined by applying the con- 
dition that ¢ must approach the correct value of 
the energy at distances far from the nuclei. 


TABLE I. Final values of coefficients relative to ¢;. 











No. of 6-point 9-point 15-point 
variations calculation calculation calculation 

a 1 1 1 
C2 0.448 0.369 0.308 
C3 .110 .100 .107 
Cs — .463 — 408 — 366 
Cs .0139 .0143 
Cs — .0082 — .0075 
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From Eq. (10) it is found that for F a finite series 


lime =—2z2?. 

TA, TBO 
Equating this to the known value of E= — 1.1026 
atomic units gives z=0.74, which is the value 
used. A slightly different value of z would not 
be expected to decrease the accuracy of the 
method. 

The coordinates of the points used are as fol- 
lows: (ra, rp) = (0.4, 1.6), (0.4, 2.4), (0.6, 1.8), 
(0.6, 2.2), (1, 1), and (1, 2.8) for the first calcu- 
lation. These six plus the following three were 
used in the second calculation: (1.4, 1.4), (1.5, 2) 
and (1, 3). For the last calculation the following 
six points were taken in addition to the above 
nine: (1.2, 1.2), (1, 2), (2, 2), (2, 3), (2, 3.5), and 
(2, 4). It should be noted that no points need be 
taken where rg <r. since the functions used are 
symmetrical in r, and rz. 

By using Eq. (10) together with (7) the 
matrices Dm», Hm,n, and Km,» were calculated 
for each of the three sets of points to an accuracy 
of eight significant figures. For the first calcula- 
tion with six points only the first four functions 
of (11) were used and the initial values of the 
coefficients were arbitrarily taken as c,=1, 
€.=0.5, c;=c,=0. Starting with c; and going on 
tO C4, C1, C2, C3, C4, etc., the coefficients were varied 
according to the procedure given in the appendix. 
After several cycles of single coefficient varia- 
tions the values relative to c, arbitrarily taken 
as 1 were as indicated in the first column in 
Table I. The gradual decrease in the mean square 
deviation with successive variations in the coef- 
ficients is shown in Fig. 1 and the successive 
changes in the value of e, in Fig. 2. For the 
9-point and 15-point calculations the variation 
procedure was initiated using approximate values 
of the coefficients obtained in the previous cal- 
culations. 

As was expected in each case as the value of 
é approaches to the correct value the mean 
square deviation of e is minimized. The actual 
deviation of about 0.01 atomic unit for the 
15-point calculation is somewhat less than would 
be expected from the lowest value of the root 
mean square deviation of 0.033. The fact that 
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the 15-point calculation does not give the best 
answer for the energy is probably due to the 
particular choice of points. There is no doubt 
that with increasing number of points and func- 
tions the correct result could be obtained as 
closely as desired. 

Application of the method to the H; problem 
is in progress. 


APPENDIX 


The Minimization of (e)s,—«* by Successive 
Approximation 


Given a set of approximate c,,’s. Calculate the 
quantities 


D= a Le CubaD an, n 


Ha) «10 Catellas (12) 
K= = i. Cobbalene 
and 
D,, = 2 a CmD m, k 
A= Yom Cn(Him, rc +A: m) (13) 


K,=2 > Calha.t 
éw=H/D, (e)w=K/D. 


Primes will be used to indicate new values of 
each quantity. Let 


C=C, + Acy.. 
Then 
D! =D+D,Ac,. + Dy.(Acy)? 


HH’ =H+HAc.+ Ayx(Ac,)? 
K'=K+K,Ac.+ Ky.(AC;)? 


(14) 


and, of course, 


é'w=H'/D’ and (e*)'y=K’/D’. 

The best value of Ac, and therefore c’, can be 
then obtained by passing a parabola through the 
resulting values of ((e*)s,— x)’ or with sufficient 
accuracy merely by inspection. The correspond- 
ing D’, H’, K’ are calculated, and using the new 
value of c,, everything is ready for the variation 
of the next C».* 


* Note added in proof: The present method is closely re- 
lated to that of D. H. Weinstein, Proc. Nat. Acad. Sci. 20, 
529 (1934), a fact called to my attention by Dr. L. Pauling. 
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Surface Tension of Micelle-Forming Solutions 


Hans M. CASSEL 
Chicago, Illinois 
February 16, 1942 


HE capillary curves of some soaps' and other long 
chain molecules in aqueous solution exhibit a peculiar 
minimum of the surface tension (s.t.) at relatively low 
concentrations (less than 0.01 molar), the s.t. becoming 
smaller than half that of water. McBain in numerous 
publications,? as well as Nutting, Long, and Harkins,’ 
believes that the observations correspond to true equilibria 
and that, according to Gibbs’ equation, the adsorbed 
quantities vanish at the concentrations of the minima. 
In other words, it is assumed that there the composition 
of the surface layers equals that of the solutions. On the 
other hand, it can be proved by surface viscosity measure- 
ments and direct tests‘ that rather dense surface films do 
exist at the same concentrations. To reconcile the con- 
tradictory conclusions, McBain has recurred to postulating 
a void of solute in the deeper layers of the solutions. On 
this basis, the existence of very thin, double-faced struc- 
tures, like soap bubbles, is hard to understand. 

In writing Gibbs’ equation: do/dm=—I-du/dm (¢ 
designates the s.t.; m, the concentration of solute; I’, its 
superficial density; u, its thermodynamic potential) one 
notices immediately that do/dm will vanish when either 
one of the factors of the second member of the equation 
tends towards zero. Since a decrease of T is out of the 
question the alternative has to be considered that a gradual 
fading of the change in potential goes on as the concentra- 
tion increases. 

Unfortunately, there is not much information available 
about the activities of these solutions in the low concentra- 
tion range and at the temperature of the s.t. measure- 
ments. Values derived by Randall, McBain, and White® 
from vapor pressure measurements show that Na palmitate 
and stearate start to approach constant activities above a 
molality of 0.3 at 90°C. Hence, it may be concluded that 
these and similar substances behave identically at room 
temperature in a much lower concentration range. 

Such a behavior is not surprising in view of the fact 
that all solutions with which the s.t. minimum is en- 
countered display definite signs, and are indeed typical 
examples, of micelle formation. Some of these solutions 
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are notoriously very close to being saturated. The condi- 
tions obtaining are analogous to certain features of mono- 
molecular films spread on water. Langmuir’ has put for- 
ward the theory that the almost horizontal portion follow- 
ing a kink in the A-F curves is due to the condensation of 
single molecules into two-dimensional micelles. 

Also J. E. Mayer’s ideas* on the condensation process 
in general seem to support the expectation that in certain 
cluster-forming systems at sufficiently low temperature the 
Gibbs potential tends to a limiting value as saturation is 
approached (Case “a” of Born and Fuchs®). This effect 
may be enhanced if the interfacial tension of the particles 
is very small.!° 

The above considerations, of course, do not account for 
the ascending branch of the capillary curve. It is probable 
that the surface films at higher concentrations assume a 
crystalline or liquid-crystalline structure. Unless such films 
are formed in a reversible way they are likely to be in a 
state of stress which, to be released, takes much longer 
periods of time than usually spent for the observations. 
A superposition of s.t. and surface stress may be reponsible 
for time effects also before the minimum concentration 
is reached. 


1H. M. Cassel, J. Am. Chem. Soc. 57, 2009 (1935). 
2 J. W. McBain and G. F. Mills, Report on Progress in Physics (1939) 


Vol. 5, & 30. 

3G. C, Nutting, F. A. Long, and W. D. Harkins, J. Am. Chem. Soc. 
62, 1496 (1940). 

4J. W. McBain and L. A. Word, Proc. Roy. Soc. A174, 286 (1940). 

5M. Randall, J. W. McBain, and A. M. White, J. Am. Chem. Soc. 
48, 2517 (1926). 

6H. V. Tartar, V. Sivertz, and R. E. Reitmeier, J. Am. Chem. Soc. 
62, 2375 (1940). 

71. Langmuir, J. Chem. Phys. 1, 756 (1933). 

8 J. E. Mayer and S. F. Harrison, J. Chem. Phys. 6, 91 (1938). 

® M. Born and K. Fuchs, Proc. Roy. Soc. A166, 391 (1938). 

10 M. Volmer in cooperation with the author, Zeits. f. physik. Chemie 
125, 156 (1927). See also J. Frenkel, J. Chem. Phys. 7, 538 (1939). 





Polarization of Luminescence in Crystals 


B. V. THOSAR 
College of Science, Nagpur, India 
February 24, 1942 


HE state of polarization! of the sharp red luminescent 
doublet of ruby at \6936, and also of the anti-Stokes 
bands shifted from the doublet by 756 cm=, 523 cm=, and 
368 cm=!, respectively, has been quantitatively investi- 
gated. The technique usually employed for studying the 
polarization of Raman lines was used, the observations 
being made visually. The crystal was kept with its optic 
axis vertical, along OZ, and the intensity of luminescence 
polarized parallel to this direction was compared with that 
polarized perpendicular to it (along OX). Table I gives the 











TABLE I. 
Lines and bands Intensity p=Iozillox 
6936 (doublet) sharp, 0.16 
intense 
Av 
756 cm7! intense, 0.70 
broad 
§23 cm=! v. intense, 0.37 
broad 
368 cm=! weak >1.00 
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ratio Ioz/Iox of the intensities polarized parallel and 
perpendicular to the optic axis of the crystal. 

When the crystal is mounted with the optic axis along 
the direction of observation OY the ratio p is unity for all 
the lines and bands listed above. 

(1) The bands have been explained as due to the super- 
position of internal vibration frequencies in the crystal on 
the sharp line radiation. Since the latter is strongly linearly 
polarized perpendicular to the crystal axis, this is a case 
analogous to the excitation of Raman lines by initially 
polarized light. The degree of polarization of the bands 
would be of help in assigning them to the different possible 
modes of vibrations within the crystal. 

(2) Regarding the polarization of the sharp doublet, it 
is to be noted that there is a striking resemblance between 
the above observations and the phenomena of polarization 
of resonance radiation of mercury, observed by Hanle.? 
The red luminescence lines of ruby behave as though the 
crystal is a natural resonance lamp with an inherent mag- 
netic field along the optic axis. The examination of the 
x-ray structure of corundum shows that the fluorescent 
Cr*** ion, which replaces the Al atom, is situated in 
crystalline fields of axial symmetry and in accordance with 
the theory* of paramagnetic crystals, the orbital moment of 
the ion should be quenched in all directions but one—the 
symmetry axis, i.e., the optic axis. The orbital moment 
being conserved only along the direction of the optic axis, 
the electric vector in the emitted radiation is confined to 
the orbital plane, i.e., the plane normal to the optic axis. 

(3) The present evidence seems to show that the Cr+++ 
ions in ruby are subjected to external crystalline fields 
strong enough to break the L-S coupling but not strong 
enough to break the coupling of individual /’s to give the 
resultant Z and of individual spins to give S. The energy 
levels of Cr*+** ion available from Bowen’s work‘ then ex- 
plain the principal emission and absorption lines according 
to the following plan: 


Transition Bowen's value Observed lines 
‘P—4F 14200 cm™ 6936 doublet (emission) 
°G—41F 15200 cm= 46580 doublet (absorption) 
*H <r} 21200 cm=! { 4752, 44762 (absorption) 
°"Da—4F 20600 cm~! \ 


This scheme is supported by the fact that ‘P—‘F transition 
being between two quartet levels is much more probable 
and accounts for the large intensity of the red emission 
lines. The other lines caused by transitions between doublet 
and quartet levels are shown only in absorption. 

(4) The occurrence of one or more lines for the *P—>‘F 
transition in Cr++* ion, when it is bedded in different 
crystals’ seems to be related either to the degree of sym- 
metry of the crystalline field in which it is situated or to 
the number of non-equivalent positions in which it can 
replace the metal atom of the ground-lattice. In spinel- 
type crystals, which are regular, there is one line and in 
AlO; and Ga2Os, there are two. This point cannot further 
be discussed here. 

I regard the above explanation of the polarization of 
luminescence of ruby as only a particular case of a more 
general interpretation of crystal luminescence, based on 
the following concepts: 
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(a) The polarization of luminescence is caused by the 
quenching of orbital moment of the “center’’ in all but 
one direction (or a few discreet ones) determined by the 
nature and symmetry of the crystalline fields. 

(b) The absence or existence of luminescence is condi- 
tioned by the complete quenching of the orbital moment 
of a possible center or its conservation along one or more 
directions within the crystal. 

Such an interpretation is being attempted. 

1B. V. Thosar, Phys. Rev. 60, 616 (1941); Phil. Mag. 26, 878 (1938) 

2 Mitchell and Zemansky, Resonance Radiation and Excited Atoms 
(Cambridge, 1934), p. 262. 

3J. H. Van Vleck, Theory of Electric and Magnetic Susceptibilities 
(Oxford, 1932), p. 288. 

4 . 


H. Van Vleck, J. Chem. Phys. 8, 793 (1940). 
5 Deutschbein, Physik. Zeits. 33, 874 (1932). 





The Migration of Colloids Through the Interface 
Between Immiscible Liquids under the 
Influence of an Electric Field, and 
the Formation of Ultra-Colloids 


G. E. GiBson AND L. BALDWIN 
University of California, Berkeley, California 
March 16, 1942 


T has been shown by Lash Miller and R. H. McPherson 

that colloidal arsenic sulfide will distribute itself be- 
tween the two phases of a mixture of water, acetone, and 
ether provided the mixture is very close to the critical 
mixing composition (plait point). A very slight excess of 
ether causes the colloid to go completely into the water- 
rich phase, when the volumetric ratios are: water : ace- 
tone : ether=1.2 : 1.9 : 1.0. 

Placing the two-phase system in a U-tube with elec- 
trodes in the upper ether-rich phases we have found that 
an e.m.f. of 200 volts will cause the arsenic sulfide to go 
back through the interface into the ether-rich phase, 
provided the ether is only slightly in excess; otherwise 
greater voltages are necessary. In similar fashion gold, 
bentonite clay, and carbon colloids have been made to 
traverse such an interface. 

The phenomenon has been viewed through the ultra- 
microscope, using gold and bentonite clay. Surprisingly 
violent motions of the colloid, not visible to the naked eye, 
became apparent soon after the application of the electric 
field. Surges back and forth parallel to the interface were 
observed, always in opposite directions in the two phases, 
accompanied by vertical wave motions near the interface. 
After some time a dark band appeared in the ether-rich 
upper phase immediately above the interface. From time 
to time particles of greater than average size approaching 
the interface in an oblique direction would shoot with 
great velocity through this dark band into the ether-rich 
phase. This dark band must be a region of high potential 
gradient, which becomes visible owing to the formation 
of an extremely fine colloid, not resolvable by the ultra- 
microscope. The formation of this ultra-colloid has been 
observed in gold and bentonite and presumably is due to 
attrition of the larger particles as they shoot through the 
interface. It probably will occur with other colloids also. 
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The violence of the surging motion is enhanced by adding 
small amounts of electrolyte (HCI to bentonite) and much 
diminished by careful dialyzing to remove electrolytes 
before use. 

We wish to acknowledge the assistance of Messrs. N. 
L. Ellis and E. Goldberg who took part in the preliminary 
experiments. 





The Absorption Spectrum of Chlorine Fluoride 


Austin L. WAHRHAFTIG 


Gates and Crellin Laboratories of Chemistry, California Institute of 
Technology, Pasadena, California 


March 16, 1942 


HE absorption spectrum of chlorine fluoride has been 

examined over the region \\3500—7000 and a band 
system found at 4800 has been studied under high dis- 
persion. These bands have single P and R branches and 
the system appears to be similar in every respect to the 
3]I9t<—!Z+ system found in iodine chloride? and iodine 
bromide.* The Deslandres Table I for the band system is 
given below. The values listed are frequencies in wave 
numbers for band origins, except those designated by h, 
which are for band heads. The bands with v’=11, 12, and 
13 are diffuse, probably because of a perturbing O* re- 
pulsive state, as in the case of iodine chloride and iodine 
bromide. 

The convergence limit for the progression v’<—0 is 
2150842 cm. The probable dissociation products are 
Cl(2P 3/2) and F (?P 1/2). The molecular constants determined 
by analysis of the bands whose origins are listed are: 
Ground state, 


D"9=2110142 cm™. 

or 60.31 kcal./mole. 
B” .=0.518+0.001 cm™. 
a’’.=0.006+0.002 cm. 
wo”. =793.2 cm, 


D” .=21495+3 cm. 


r”’-=1.625A. 
r’’9=1.630A. 
wx” =9.9 cm. 


Excited state, 
D')=2790+50 cm™. 
B’.=0.372 cm". 

a’. =0.014 cm™. 

r’-=1.92A. 

w’e= 313.484 cm. 
w’ ¢y’e= —0.400 cm. 


D", and D”’, were calculated assuming the dissociation 
products to be as given above. B’., a’, and the excited 
state vibrational constants represent only the data for 
v<1i1; that is, below the perturbation region. 


D’.=2946+50 cm“. 
B’ 14=0.132 cm", 
B’;;=0.112 cm}, 


w’ x’ =2.217 cm. 
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TABLE I. Deslandres table for the band system at 44800 











v v’ =0 v’=1 
6 20395.82 
7 20619.32 
8 20820.37 20047 .04 
9 20996.94 
10 21145.54 20372.07 
11 h 21255.5 
12 h 21338.3 
13 h 21398.4 
14 21445.87 
15 21479.95 
16 h 21499.8 
17 hk 21507.8 
Isotope bands 
8 20808.6 
9 = 
10 21134.3 








The dissociation energy of chlorine‘ is 57.20 kcal./mole. 
From this and the dissociation energy of chlorine fluoride 
listed above, one obtains: 


Cl+2F->2CIF AH=—63.42 kcal. 
or 
Do(F2) +2Q,(CIF) = 63.42 kcal., 


in which Do(F2) is the dissociation energy of fluorine and 
Q;(CIF) is the heat of formation of chlorine fluoride, 
rigorously at O°K, but with very little error at room 
temperature. 

The values given in the literature for Do(F2), all de- 
termined indirectly, range from 62 to 65  kcal./mole.’ 
Ruff and Menzel® report Q;(CIF)=27.4 kcal.; Freden- 
hagen and Krefft? state that chlorine and fluorine form 
an explosive mixture, and that the reaction is more exo- 
thermic than that between chlorine and hydrogen, for 
which Q;(HCl) = 22 kcal. It is apparent from the equation 
above that at least one of these values is greatly in error. 
At present it seems most reasonable to assume that 
Q;(CIF) is much smaller than 27 kcal., but appreciably 
greater than zero, and that Do(F2) is somewhat smaller than 
has been generally believed. 

A complete report on this research will be published 
when time permits. 

1W. G. Brown and G. E. Gibson, Phys. Rev. 40, 529 (1932). 

2J. H. Van Vleck, Phys. Rev. 40, 544 (1932). 

3 W. G. Brown, Phys. Rev. 42, 355 (1932). 

4G. Herzberg, Molecular Spectra and Molecular Structure (Prentice- 
Hall, 1939), p. 485. 


5 See, for example, reference 4, p. 486, and W. M. Latimer and J. H. 
Hildebrand, Reference Book of Inorganic Chemistry (Macmillan, 1940). 


p. 154. 
6 Otto Ruff and Walter Menzel, Zeits. f. anorg. allgem. Chemie 198, 
375 (1931). 

7K. Fredenhagen and O. T. Keefft, Zeits. f. physik. Chemie A141, 
221 (1929). 
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